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ABSTRACT 



POLARIZATION OF THE COSMIC BLACKBODY RADIATION 

AT 3.2 CM. 



GEORGE PETER NANOS, JR. 

* 

A 3.2 Dicke Radiometer configured as a polarimeter was 
used to make measurements of the linear polarization of the 
Primeval Fireball along a declination of 40.35° N. The 
resolution of the instrument was 15° in declination and 1 hour 
in right ascension. Harmonic analysis of the Stokes parameters 
of the radiation, Q and U, sets limits on sidereal variations, 
which in turn limit the possible contributions of cosmological 
origin. Results are related to the isotropy of the cosmic 
background radiation and the symmetry of the universe. 
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I. INTRODUCTION 



1.1 HISTORY OF THE FIREBALL 

In I 929 it became clear from the work of Edwin Hubble that the 

universe is expanding , and that the recession velocity of a galaxy 

is nearly a linear function of its distance.^ This result gave 

impetus to the development of expanding universe cosmologies based 

on solutions to Einstein T s equations 0 Of the many ensuing models, 

among the most successful and well established is the homogeneous 

and isotropic universe with cosmological constant equal to zero, 

2 3 

based on the work of Friedman and Lemaitre. 5 

One important feature of this model is the existence of a 
singularity at t=0; i.e., at some finite time in the past the scale 
factor in the universe was small and the density of matter and 
radiation was very high. This feature is known as the "fireball" 
hypothesis and the subsequent expansion and evolution of the 
universe is called the "Big Bang". 

During the early stages of the fireball expansion, the matter 
in the universe was extremely hot, much greater than 10^ °K, and 
was in thermal equilibrium with the radiation. Being thermal, this 
radiation had a Planckian spectrum given by 

8nhv 1 ^ 

dE. = 3 hu/ kT dv l* 

1 c ell 

The first description of this type of blackbody-radiation filled 
model v/as given by Tolman in 1931 1 Later he shov?ed that under 
adiabatic expansion the radiation spectrum is still given by (l.l), 



2 



where the tempera ture T is now given by 



T ' (t) = T 



o a(t) ^o 1+z ’ 



( 1 . 2 ) 



where a(t) is the scale factor of the universe at time t, and 1+z 

is the redshift due to the first order Doppler effect from the 

expansion, viewed by a comoving observer. 

Alpher and Gamow in 19^8 needed the fireball in their attempt 

to explain element formation in the early universe and in their theory 

derived a temperature of 5 °K for the remnants observed today, a figure 

o 5 

quite close to the 2.7 K obtained experimentally. Interest in this 
result waned as other explanations of element formation proved more 
attractive until, working independently, R. H. Dicke and his colleagues 
at Princeton were drawn to the idea to explain the reprocessing of 

g 

heavy elements into hydrogen in a closed oscillating universe. 

At the same time that two of these co-workers, P. G. Roll and 
D. T. Wilkinson, were building a radiometer to search for this 
radiation, Penzias and Wilson at Bell Laboratories realized that 
what they previously had considered to be excess receiver noise at 

7.35 cm. and had tried to eliminate, was in fact the remnant of the 

7 

primeval fireball that the Princeton group was searching for. 

Soon afterward the 3.2 cm. measurement of Roll and Wilkinson con- 
firmed the result and showed that the spectrum over this range was 

o 8 

consistent with that of a 2.7 °K blackbody. Following this, a host 
of other measurements have been performed^ ^ all of which tend to 
support the hypothesis to a greater or lesser degree, at least in 
the Rayleigh- Jeans part of the spectrum. At wavelengths greater than 
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a few tens of centimeters or shorter than a millimeter, the situation is 

not as clear. At longer wavelengths contributions due to sources 

within the galaxy make separating out the blackbody difficult while 

short wavelength observations suffer greatly from absorption and 

emission in the atmosphere. If one excludes the rocket data at very 

short wavelengths , the reliability of which seems open to question 

at the present time, Peebles has shown that jthe available data supports 

29 

the fireball reasonably well. 

ANISOTROPY 

Accepting the idea that one is looking at fossil radiation from 
early in the history of the universe, it is natural to ask what 
information it carried from this early epoch. As the early universe 
expanded and cooled, at one point the primeval plasma recombined and 
the mean free path for photon scattering became very large, so 
large in fact that most of the fireball photons we observe were 
probably last scattered during the recombination at z=1000. Only 
a subsequent reionization of the matter could bring the scattering 
epoch closer. Therefore, the spatial distribution reflects the 
distribution and Lemperature of matter at this early time. 

Tests for isotropy in the distribution of the radiation have 
been conducted by several investigators including the Princeton 

_ oQ 

group. At the present time no significant departure from 

spherical symmetry has been found for angular resolutions on the 
order of 1° of arc. In addition to early symmetry of the fireball, 
the data can be used to deduce our peculiar velocity relative to the 
local comoving frame. This effect appears as a 2h hr. asymmetry with 

39 -hi 

a hot spot in the direction of motion due to the added Doppler shift. 
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1.2 COSMOLOGICAL CAUSES OF POLARIZATION 

42 

As was pointed out by Martin Rees in 1968, the type of expansion 
asymmetries which would give rise to anisotropies in the microwave back- 
ground would also give rise to linear polarization of the radiation 
through the mechanism of Thomson scattering. Following his lead 
and using much the same notation, I will briefly demonstrate the 
existence of a linearly polarized component pf the background, given 
an asymmetric expansion. 

In this I will consider only one particular class of asymmetric 
universes, those with axially symmetric Euclidean metrics of the 
form ds^ = dt^ - A(t)^ (dx^ + dy^) - W ( t ) ^ dz^ and containing no 
magnetic field. Models of this type have been investigated extensively 
by Thorne (1967),^ Jacobs (1968),^ and others and are particularly 
attractive because of their simplicity, their close relation to the 
standard Einstein-de Sitter universe, and the straight forward way 
in which they lead to a polarization. In the case of three 
inequivalent directions, the situation would be much more complex 
(there being azimuthal dependence to the equations) and the calcula- 
tions correspondingly more difficult. 

Taking our given metric we can first define Hubble's constants in 
each direction, a = dA/Adt and w - dW/Wdt. From these we can compute 
an average expansion constant h = (2a+w)/3 and a differential constant 
Ah = (w-a) and using them, can discuss the existence of a temperature 
asymmetry in this type of universe. As was mentioned in the preceding 
section, Tolman showed that in the Einstein-de Sitter case the 
temperature of radiation is inversely proportional to the scale 
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factor, which in turn is proportional to V^"/ 3 . In our case V the 
volume is proportional to aS/ and in a similar way we can write 

T a (A 1 /3 j (1.3) 

where T is an average temperature over all directions of observation. 

At any time t, one can define an approximate isotropy parameter equal 
to AT/T = e. AT is the temperature difference between the A and W 
directions and can be written ? 



AT = t d/dt (T - T ) 
s w a 

a t d/dt (l/W - l/A) 
s 

~ t Ah/(A 2 \-/) 1//3 . 

— s 

Here t is the time since the last scattering, and is given by a , 
s 

the mean free path of the photons , divided by the velocity of light ; 

i .e. 



t 

s 




1 

nl(t)o T c 



(1.5) 



where n = the number density of electrons 

l(t) = the percentage of them which are ionized 
= the Thomson scattering cross section. 

Putting this form into the previous equation, we have finally, 



e = (A^V/f ' 3 A T 
Ah 

*I(t)o T c ‘ 

To obtain a more precise result requires consideration of the 
cosmological redshift in an asymmetrically expanding universe as is 
done in Appendix A2. In particular for the axial model being 
considered, the temperature of the radiation as a function of the 
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direction of observation with respect to the axis of the universe is 
shown to be 

A ^ W 2 i 

T = T [-f- sin 2 0 + cos 2 0]" 2 , (1.7) 

A W 

in the absence of scattering. 

Throughout this discussion the zero subscript will refer to 
the value of a quantity at the present time.^ In the event that the 
expansion asymmetry is small, as it is in the case we are considering, 
this equation can be cast in a little more convenient form as 

T - T (1 + e sin 2 9). (1.8) 

Taking this type of temperature distribution as a starting point, 
the polarization of the microwave background through the mechanism of 
Thomson scattering can now be demonstrated. As shown in Figure 1.1, 
we will define the polarization of an electromagnetic wave from the 
microwave background as seen by an observer in terms of the plane 
containing the line of sight to the observer and the direction of 
asymmetric expansion. The direction perpendicular to this plane is 
called the 'a' direction while the direction which lies in this plane 
and is, of course, perpendicular to the direction of propagation is 
called the ccplanar or 'w' direction. Following Rees, I will 
write the temperature of the radiation polarized in each direction in 
the following way: 

T = T (1 + e sin 2 0) 

T w = T (1 + e w sin 2 0) , (1-9) 

where T now refers to T 

av 
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FIG. \ A 



DEFINITION OF THE POLARIZATION REFERENCE 
PLANE 
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By doing so, I can characterize the polarization completely by means of 
the parameters e and e . In fact, the evolution of the polarization 
with time may be ascertained by constructing and solving the differen- 
tial equation, describing the change of these parameters with expansion 
and scattering of radiation. It must be understood that this type of 

approach is only valid when the temperature is low enough, say 
Q o 

10 for elastic Thomson scattering to he the primary mode of 

interaction with the photons and when the asymmetry in the expansion 
is small enough that the above expansions of the temperature distribu- 
tion are valid. In light of the small observed temperature asymmetry 
this latter condition doesn't seem too unreasonable. At this time there 

doesn't seem to be any evidence for a very hot plasma greater than 
9 o 

T = 10 K at a recent epoch . 

If we express the temperature of the radiation in the two 
polarizations of the microwave background in the way described above, 
we can compute the effect on the asymmetry parameters e & and of 

a single Thompson scattering. In Appendix A3 the result in matrix 
form after integration over all angles of incidence -is shown to be 



C-) • (. 

w 



0 0 ' 
1/2 7/30 > 



C) ■ 

W 



( 1 . 10 ) 



where the primed quantities refer to the values after scattering. 

The most notable effect in this calculation is that the asymmetry 
in the orthogonal 'a' component is destroyed by the scattering process 
while the asymmetry parameter in the coplanar polarization is preserved. 
If to begin with, both parameters had been equal to e , the effect 
of the scattering would have been to increase the maximum polarization 
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from zero to l/30 e, which demonstrates the effect of scattering in 
producing a polarization. 

Because we know the time scale over which this change occurs, 
we can write an expression for the rate of change of the asymmetry 
parameters with time as follow's : 



d_ 

dt 









( 1 . 11 ) 



In addition to the effect of scattering on e and e , one must 

a w 

consider the effect of expansion even though this is not what leads 
to the existence of a polarization. From equations 1.5 and 1.6 we 
find that the growth of a temperature asymmetry with time due solely 
to expansion is t g Ah. This implies 



d_ 

dt 





At 



t s Al1 jj) 
t 



= Ah (\) . 



( 1 . 12 ) 



s 

Combining this with the previous result w r e obtain a final differential 

equation for e and e : 

^ a w 



ft ( e w) = ^ (i 



) + 



nl(t)o T c 



-1 0 
-1/2 -23/30 



%)■ (1.13) 



Using the result from Appendix A1 that Ah =Ah^/A^W and noting that 
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n = n. 



0 



/A^\l simplifies the expression to 



d_ 

dt 




n Q l(t)a T c 

aSj 




-1 0 
-1/2 -23/30 




• (1.14) 



It is also shown in Appendix A1 that A and W can be expressed in 
the following way : 



which demonstrates their time dependence explicitly, leaving only one 
unknown function of time in the equation l(t). The use of q here 
should not be confused with the deceleration parameter. I only use 
it for consistency with Thorne's notation. Since we have no 
theoretical or experimental handle on the functional dependence of 
I(z), I will follow Rees and set it equal to one in order to obtain a 
useful approximate solution. Adding this to the above, we obtain the 
following two differential equations for the asymmetry parameters: 

de n c m c Ah 

o r P rs 



A - A (t/t ) 2 / 3 
o w o' 

W = A + qA -1 / 2 , 



(1.15) 



a 



dt 





de 

w 




Ah 



o 



dt + .3 



A 3 (t/t ) 2 +qA 3 / 2 (t/t ). 
o ' o o ' ' o 



(1.16) 



Solutions to these are easy to obtain and are given by 
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All 






6 a n ct_c 
o T 



1 + 



n a ni ct /qA 
o T o' ^ o 



3/2 



15Ali 



A 3 /2 t 

O / 



at 



e = 
w 



w 



3/2 



n a m ct / A 
o l o T o/q o 

1 + ~ 3/2 t I 



(1.17) 



which asymptotically approach 
Ah 



e = 

a 

11 °m C 

o T 
15Ah 

e = o 
w ~ 

7n o°T° 



(1.18) 



In the event that q is a small number, as long as t > qt^/A 
these values are reached and, in fact, it is shown in Appendix A2 
that e being small is equivalent to requiring that a be small. 

From the asymptotic solutions we can now compute e = \ (e a + e^) 
and the degree of polarization (T - T^)/T. 



e = 



ni Ah 
11 o 

7 n o a T o 



= e 



(T - T )/T = 
v w a y/ 



8 

7n o CT T C 

- 8 - I 
11 



(1.19) 



2 

sin 0 



• 2 „ 
sm 0 



This is the scattering induced polarization which we wished to 
demonstrate. Although the assumption that l(t) is identically equal to 
one is unrealistic, the solution so obtained may now be applied to more 
tenable pictures of the evolution of the universe to test for an 
anisotropic expansion. 
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II. APPARATUS 

2.1 THE MEASUREMENT OF LINEARLY POLARIZED RADIATION 

In order to measure and understand the effect of a cosmological 
linear polarization, we must first understand ho >w polarized radiation 
is measured. The simplest kind of polarimeter is shown in Figure 2.1a. 
It consists of a polarization sensitive anterma, of which the simplest 
is the center fed, linear dipole shown, followed by a radio receiver. 

To measure polarization, you aim the polarimeter at a source and rotate 
it until you obtain the maximum signal. The angle between the axis of 
the dipole and some suitable reference position is the angle of the 
polarized signal; and the difference between this signal and lowest 
signal received, as the device is rotated through 360 degrees, is 
proportional to the magnitude of the polarization. This method is 
fine and, because of its simplicity, is probably the best as long as 
the signal is strong enough to give a clear difference between maximum 
and minimum. 

In many practical applications the signal is quite small and has 
to compete with a large amount of system noise from the receiving 
equipment itself. In this case the apparatus in Figure 2.1b is 
better. It consists of two antennas of the dipole type with their 
polarization sensitive directions oriented at right angles to each 
other and each connected to its o wn receiver. When this is aimed at 
a source, a signal S_^ will appear at the output of each receiver: 

5 1 * c < T u 4 T p c °= 2e * T RECl > 

5 2 = c < T u * T P si " 2e 4 T r,EC ? >. 



( 2 . 1 ) 
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where is the temperature of the unpolarized radiation from the source, 
T is the temperature of the polarized component, 

6 is the angle of the polarization with respect to antenna one, 

C is a proportionality constant depending on the gain of the 

receivers, the method of detection, ' and the receiver band- 
widths . 

T is the equi valent noise temperature of each receiver. 

REC * 

Subtracting the two signals gives an expression for and 0 3 



V S L' = c[T p cos2e + < T BEC - T REC„ )] - 



( 2 . 2 .) 



'1 * — 2 

If the noise figures of the two receivers are the same, the second 
expression in the brackets has mean zero and, for sampling times long 
compared to the reciprocal of the receiver bandwidth, a Gaussian 
distribution. Thus the time average of this term is zero and we can 
define 

< S, - S 0 > m ._ 

(2.3) 



/\ < ^1 °2 ^ime m „ on 

Q = = T cos20. 

c p 



Rotating the antenna system by forty- five degrees gives another 
pair of signals. 



S l’ = C ^V T P COS ( 0 + n /4) + t rec ^ 

V = c[ V T P sin 2 (0 ♦ n/4) ♦ T BEC2 ], 

from which a second parameter can be defined in a similar way 



(2.4) 



U = < sp- S ? ,’ > = T sinaa. 



(2.5) 



In terms of these parameters the polarization temperature and angle 



are given as 
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ANTENNA 



METER 



(a) 




(b) 

FIG. 2.1 



TWO TYPES OF POLARIMETER 



( 2 . 6 ) 



0=1 ARC TAN (|) . 

It is clear that these two parameters, known as the Stokes 
parameters, uniquely determine the polarization • of a source and thus 
together form a complete set of orthogonal measurements to determine 
this quantity. 

* 

For many applications this type of polarimeter works quite well, 

but it too has limitations. Even when one finds two receivers with 

well matched noise figures, drifts in the gain of the devices, which 

are not necessarily random in nature, can cause problems. If the 

signal one is seeking to measure is of the same order as the product of 

the receiver noise and the size of the average fluctuations in gain 

between the two receivers, serious errors in the measurement result. 

These errors cannot be averaged out. A solution to this problem, 

diagrammed in Figure 2.2, is the replacement of the two receivers 

with a single Dicke receiver which is set to switch synchronously 

45 

between the two antennas . The power entering each antenna is 
written as 



= T + 


2 

T cos 9 




u 


P 

2 


(2.7) 


= T + 


T sin 9. 




u 


P 





The signal after the switch and before the receiver is given by 

Q = T x - T 2 = T p cos20 f(uit), (2.8) 

where f(urt) is a periodic function at the switching frequency, with 
amplitude one, v/hich is dependent upon the switching characteristics of 
the receiver and the modulating signal. After passing through the 



l6 




h" 



C\i 

OJ 

d 

u_ 



< 



DICKE POLARIMETER 
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receiver, the signal picks up a contribution due to receiver noise, which 
is a sum of contributions at all frequencies up to the cutoff of the IF 
and video amplifiers. In parti cular, the contributions due to gain 
fluctuations are at low frequencies , which means that by choosing the 
switching frequency high enough, we can isolate the wanted signal from 
the gain fluctuations. Now if the signal is processed through a phase 
sensitive detector like a lock-in amplifier , only the wanted signal 
and those noise contributions with frequencies in a limited bandv/idtli 
around the switching frequency are recovered and the contributions 
from gain fluctuations are excluded: 

S = c[Q + (Av)]. (2.9) 



Again taking a time average we can recover the parameter Q, 



<C> 

Q = Time . (2.10) 

c 

If the antennas are rotated by 45 degrees, U can be measured in a 
similar way. 

The above progress is not gained without cost. By switching between 
the two antennas we are only looking at each one half, of the time and 
effectively throwing away half the available signal. The minimum 
detectable signal in a radio receiver as a function of system noise 

46 

temperature, integration time, and bandwidth is given by 



A T . 



rain 



T 

REG 

/AvAt 



( 2 . 11 ) 



This tells us that if we have a receiver with bandwidth Av and we wish 
to measure Q or U to an accuracy AT then we must carry out the 
above time average for a period At in order to eliminate the effects 
of For the Dicke receiver, since we are throwing away half the 
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signal, things become worse by a factor of two: 



A 



T . 
min 



T 

2 REC 
/AvAt 



( 2 . 12 ) 



If, instead of a wide band pre-multiplication amplifier, a narrow band 
tuned amplifier is used, then the first harmonic at the switching 
frequency alone is extracted and the corresponding loss in sensitivity 
is 4/tt /2 . This loss is more than made up i'or by the freedom from 

overloading and simpler amplifiers made possible by this approach. 
Finally, if the switching is also done sinusoidally (something which 
is often done as a matter of convenience), the sensitivity is further 
reduced by the factor 4/tt. Putting these together gives the expression 
for a sinusoidally switched Dicke radiometer with narrow tuned video 
amplifier 



A 



T . 
mm 



2/2 



T 

REC 



(2.13) 



/ AvAt 
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2.2 TERRESTRIAL OBSERVATION OF A COSMOLOGICAL LINEAR POLARIZATION 

Next we ask, what kind of signal does polarized background 
radiation produce in an earthbound radiometer, given the orientation 
of the axis of the asymmetric expansion with respect to the rotation 
axis of the earth, the latitude of the observer, and the direction 
in which the observer is looking. In order to simplify the question 
and the subsequent analysis, I will assume that the earthbound 
observer will make all his observations by pointing his polarimeter 
at the zenith. It will turn out later that there are sound experimental 
reasons for this procedure in addition to the computational ones. 

If we were to step away from the earth for a moment and aim our 
polarimeter in a direction at an angle 0 with respect to the axis 
of the universe, and if we aligned our polarimeter with the sensitive 
axis of the antenna aligned with the coplanar direction, then from 
equations 1.9 v/e find that we would measure a Stokes parameter 0, 
given by 



or in an equivalent, and for ini'- purposes more convenient, form as 



Nov;, stepping back on the earth which is rotating about an axis, which 
forms an angle § with respect to the symmetry axis of the universe, I 
observe the zenith which is a direction whose declination is given by 
the latitude of the point of observation and whose right ascension is 
just the local sidereal time at the observation site. Figure 2.3 is a 
schematic representation of the situation, showing the experimental 




(2.14) 
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T - T = (T - T ) sin 9. 
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SYMMETRIC COSMOLOGICAL POLARIZA- 
TION AND THE ROTATION AXIS OF 
THE EARTH 
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Given this situation, I show in Appendix A4 that the polarimeter 
will measure a quantity 

T.. - T 0 = (T - T ) )[cos^6(l - ^ sin^§) + sin^§(l - cos^6)cos2(g - t) 

1 2 v w a'max I v 2 ' v 2 ’ v o ' 

' 2 2 

- l/2 sin2? sin26 cos(q' o ~ t)] cos p - [sin § sin 6 sin2 (or - t) 



2 ( 

- sin2$ cos 6 sin(o Q - t)] sin p) , 



(2.16) 



where t is the local mean sidereal time at the observation site, 

1 

p is the angle between the sensitive direction of the antenna 
and the direction of the Worth Celestial pole, 

6 is the declination of the zenith, and 

is the right ascension of the symmetry axis of the universe. 

As the polarimeter sweeps out a circle of constant declination night 
after night, we can add the output records synchronously at the sidereal 
rate and, in this way, perform the time average needed to overcome 
the noise temperature of the receiver. If we alternate the orientation 
of the radiometer between two positions 45 degrees apart, we can make 
the two orthogonal measurements, necessary to uniquely specify the 
polarization, simultaneously. 

From equation 2.l6 we can see that by grouping all the constants 
together, the result can be expressed more simply as 

T, - T_ = (T - T ) /a + Bcos2(g - t + 6 n ) + C cog (g - t + 6 0 ) i , 

(2.17) 



with relationships between the coefficients easily derivable from the 
parent equation. This form displays the harmonic dependence upon local 
sidereal time and thus the dependence upon right ascension. In general 
then, each measured parameter will have a 12 hour and a 24 hour harmonic 
component. It is the object of this experiment to measure the polariza- 
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tion parameters as well as possible and look for the harmonic dependence 
which would indicate the existence of a large scale polarization of 
cosmological origin, obeying this simple model. 



2.3 THE FARADAY SWITCHED POLARIMETER 

The device used to perform this experiment is an offspring of the 
Dicke type receiver called a Faraday switched, polarimeter . As is 
shown in Figure 2.4, instead of two polarization sensitive antennas 
connected to the receiver through a single pole double-throw switch, 
one conical type polarization insensitive antenna is used. This horn 
feeds into an axially symmetric chamber containing a thin ferrite rod. 
Application of a sinusoidal magnetic field to the ferrite causes the 
incoming radiation to be rotated back and forth synchronously. The 
chamber is followed by a dual mode transducer which resolves incoming 
radiation into two orthogonal components and thus acts as an analyzer. 
One port of the dual mode transducer is connected to the receiver 
while the other is dumped into a suitable matched load. The net 
effect of this arrangement is the same as was obtained using two 
antennas, only now some of the matching problems have been eliminated. 
With the field on the ferrite in one direction, radiation from one 
polarization is rotated by the Faraday effect and passes through the 
dual mode transducer into the receiver. If the field is then reversed 
and the magnitude has been adjusted so that the radiation is rotated 
90 degrees in going from one state to the other, the orthogonal 
polarization will now be presented to the receiver. In this way a 
switched signal, proportional to the temperature difference between 
the two directions, is produced. 



CIRCULATOR 








FARADAY SWITCHED POLARIMETER 

FIG. 2.4 
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The rest of this chapter will be devoted to this de-vice and its 
calibration. Frequent reference will be made to Appendix B where 
the detailed specifications and measured parameters are given. 

2 . 4 THE ANTENNA 

A conical optimum gain horn with an aperture choke for suppression 
of the backlobes was used in this experiment.* Provision was made for 
insertion of a calibration probe in the side of the horn and this 
feature was used for part of the time. Later analysis showed that the 
probe increased the sensitivity of the polarimeter to ground radiation, 
an effect which is described quantitatively in Appendix Bl. 

The full width of the main lobe between half power points was 
about 15 degrees, giving the polarimeter a resolution of one hour 
in right ascension. 

2.5 THE SWITCH 

The heart of the polarimeter is the Faraday switch used to rotate 
between orthogonal polarizations. Because it was used to produce the 
signal, any departure from ideal operation will intimately affect the 
results of the experiment. In this section I will present those 
principles of operation necessary to an understanding of the polarimeter. 
A more detailed analysis of the operating principles and parameters of 
the switch and how they are affected by outside influences is given 
in Appendix B2. 

The basic principles can be most easily understood by considering 
the switch as an ideal Faraday rotator. As shown in Figure 2.5, the 
active element is a thin cylindrical piece of ferrite mounted coaxially 
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THE SWITCH 



FIG. 2.5 
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in a section of circular waveguide around which is wrapped a coil of 
wire. Putting a current through the coil produces a longitudinal 
magnetic field in the ferrite which in turn produces the Faraday 
effect, by causing the magnetic moments of the ferromagnetic ions to 
align with the field. If then, a circularly polarized wave propagates 
along the axis, its velocity is either decreased or increased depending 
on whether its helicity has the same or opposite sense, relative to 
the precession of the ferromagnetic moments. Since a plane polarized 
wave is composed of equal left and right circularly polarized components, 
one component is retarded relative to the other and the plane of 
polarization is rotated. By placing an oscillating current of 
sufficient magnitude on the driving coil, the propagating wave can be 
rotated first in one direction by 45 degrees over the length of the 
ferrite, and when the current is reversed, 45 degrees in the opposite 
direction. Now if the rotator is followed by a dual mode transducer 
which acts as an analyzer, separating the radiation into two orthogonal 
components, one obtains, in either arm of the transducer, radiation 
which is modulated at the switching frequency. The amplitude of the 
modulated signal is proportional to the difference between two 
polarization components of the incoming radiation, referenced to a 
coordinate system, which has one of its axes bisecting the angle 
between the sensitive directions of the dual mode transducer. The 
graphs in Figure 2.6 depict the result for incident polarized radiation 
in two situations. In one ease, the direction of polarization is 
parallel to one of the axes of the reference coordinate system, giving 
the maximum signal. In the other case, the incident radiation is at 
an angle of 45 degrees to this direction or, in other words, parallel 
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FIG. 2.6 

SIGNAL PRODUCED BY THE SWITCH FOR 
TWO ORIENTATIONS OF INCOMING 
RADIATION 
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to the direction of one of the sensitive directions of the dual mode 
transducer, giving zero signal. That this second case in fact gives 
zero signal comes about because the output has a periodicity at twice 
the switching frequency. Since the lock-in extracts only the first 
harmonic of the switching frequency, this contribution is zero. 

It should be noted that I am looking at temperature, a quantity 
which for an extended blackbody source is related to the incident 
power by 

P = kTAv, (2.18) 

in the Rayleigh-Jeans portion of the spectrum. Therefore, I am 
sensitive to a quantity which is proportional to the square of the 
incident field. Because of the lack of dependence on the sign of the 
incident field, positions of the polarimeter which differ by a 
rotation of l 80 degrees about the ferrite axis are equivalent, while 
positions 90 degrees apart are opposite in sign. 

In practice this symmetry is not perfect, being broken by an 
offset due to absorption and emission of radiation by the switch 
itself. To make matters worse, this offset is a function of three 
major perturbing influences; the temperature of the ferrite, changes 
in ferrite magnitization due to DC fields like the earth’s magnetic 
field, and mechanical strains transmitted to the ferrite through its 
mountings. In order to reduce the effects of temperature, the 
polarimeter was mounted in a temperature controlled environment, 
where, except under extreme conditions of heating or cooling, the 
temperature was maintained to about one degree centigrade. Effects 
from magnetic fields were removed by encasing the entire microwave 
front end of the device in a shield constructed of extremely high 



permeability material. Finally, mechanical strains were kept to a 
minimum by mounting the axis of the pclnrimeter vertically, and 
thus parallel to the gravitational force. 
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One can choose to take the signal out of either arm of the dual 
mode transducer, since these outputs give the same signal, differing 
only by 180 degrees in phase at the switching frequency. In fact, 
it would be possible to use both outputs to advantage as will be 
discussed later. In this experiment I used only one arm, attaching 
a small standard gain horn, aimed at the sky, to the other. In this 
way the unwanted signal was discarded and the unused port of the 
switch was terminated in a cold load. 

Returning to the subject of the switch offset, though I have 
described how perturbing influences were kept from altering it, I 
haven't yet explained how to eliminate it as an instrumental effect. 
The easiest way to accomplish this is to take alternate readings of 
the polarization, separated by a rotation of 90 degrees about the 
ferrite axis. This would result in two signals 



S n = T. + T __ , 

1 1 offset 



S 2 “ " T 1 + ^offset 5 



and by subtracting them we would obtain 



S l" S 2 = 2V 



(2.19) 



( 2 . 20 ) 



a quantity in which the offset has been subtracted out. Of course, 
this procedure depends on being able to shield out the effects of the 
earth's magnetic field and the perturbing influence of gravitational 
strains. Having the polarimeter axis vertical keeps the gravitational 



forces constant with rotation. 
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2.6 THE CIRCULATOR AND ISOLATOR 

These components are designed to provide isolation between the 
microw 7 ave mixer and the back of the switch. Because the local 
oscillator is not perfectly monochromatic and still has a power 
contribution at the sideband frequencies, isolation is made necessary 
to prevent power leakage back up the waveguide. For the Gunn 
oscillator used, the power in the sidebands yas on the order of 
-l40 db below the center frequency power. As this is a factory claim 
and wasn't directly measured, adding a factor of ten might be 
appropriate; so -130 db might be more comfortable. For input powers 
of between one and two millawatts this means that there may be as 
much as 2 x 10 watts left in the sidebands. It might not seem like 
much but being in the Rayleigh-Jeans portion of the spectrum and with 
a bandwidth of 96 MHZ, we find that a millidegree corresponds to 

13 

1.4 x 10 watts. In short, we may find as much as 100 millidegrees 
worth of power leaking out of the mixer in the sidebands we are 
sensitive to. 

Ordinarily this would probably be no problem, but the reflection 
coefficient of the back of the switch is a strong function of magneti- 
zation. For that reason one doesn't want to provide it with a nice 
source of power at the sideband frequencies, which it could modulate 
through reflection and send back, masquerading as real signal. 

To provide isolation, a circulator and isolator were inserted 
between the switch and mixer. Figure 2.7 diagrams the operating 
characteristics of the circulator. Power leaving the switch is passed 
to the mixer. Power from the mixer is passed to the third port which 
is terminated in a standard gain horn aimed at the sky, and a cold 
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CIRCULATOR 



FIG. 2.7 
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5 °K from the sky illuminates the back of the switch, giving it little 
signal to modulate by reflection. About 30 db worth of isolation is 
provided in this way. 

Another 40 db is provided by the isolator, which is inserted 

between the circulator and the mixer. It, in fact, is also a circulator 

but with its third port terminated internally in a matched load. The 

addition of these two devices reduces the residual local oscillator 

-23 -5 

power reaching the switch to about 10 watts or 10 mdeg. K> 
a much more acceptable level. 

Because these two devices depend on ferrites for their operation, 
it was thought that perhaps, like the switch, they would be subject 
to variations in their operation due to stray magnetic fields. It 
might not have been as severe in this case, as neither of these 
devices uses an oscillating field to produce a signal. None the less, 
they pass the offset produced by the switch and if their absorption is 
changed by a large enough amount, they could produce an orientation 
dependent signal. For this reason they were tested by an application 
of field along the axes of their ferrite devices. This should be 
their most sensitive directions. No measurable effect was observed 



for fields up to 11 Gauss. 
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2.7 THE RECEIVER 

A standard off the shelf model heterodyne receiver, having a 
separate double balance mixer and IF amplifier was used at 9*37 GHz. 
The width of each sideband was 48 MHz , giving a total receiver band- 
width of 96 MHz. Local oscillator power was provided by a Gunn 
oscillator coupled through a frequency meter and attenuation pad. 

The only unusual operating characteristic of this receiver was 
a bad impedance match between the mixer preamplifier and IF amplifier. 
If a resonant Length of cable were used and the system perturbed by 
some electrical interference, oscillatory breakdown could occur. 

By judicious choice of cable lengths and careful shielding against 
interference, the problem was eliminated. Because a chart record of 
the lock-in output was kept and because this problem gave a striking 
signature, those times that it did occur (as the result of extreme 
electrical interference in thunderstorms, for example) were easily 
seen. and the affected data eliminated. 

One modification was made to the standard IF amplifier and that 
was the addition of a BNC terminal to allow measuring the DC voltage 
at the second detector before it was capacitively coupled to the 
video stage. As will be explained in the next section, this was 
needed for the calibration procedure used. 
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2.8 CALIBRATION 

Because there were no standard polarized blackbody sources 
available for calibration at 3.2 cm, another method had to be employed. 
The procedure used involves two steps. First, the change in voltage 
at the second detector as a function of antenna temperature is deter- 
mined, followed by a measurement of the gain through the rest of 
the system. * 

The calibration was begun by measuring the DC voltage at the 
second detector with the polarimeter looking at the sky. The size of 
this voltage is proportional to the power entering the horn from the 
sky plus some offset due to receiver noise and biasing of the 
electronics. Next, a piece of eccosorb microwave absorber, a good 
unpolarized blackbody source at ambient temperatures, was placed over 
the antenna and the voltage again read. This time the voltage is 
proportional to the power from the eccosorb plus the same offsets 
that were present before. In this way we have measured the following 
two quantities : 

V = c (T + T ) 

sky 1 v sky offset ‘ 

V . = C n (T . + T „ ,), (2.21 

eccosorb 1 N eccosorb offset 7 ? v 

where is a constant dependent on the gain of the receiver. Combina- 
tion of these two allows one to determine and thus the response 
in Volts/°K of the receiver to this point. Now imagine that the 
polarimeter is looking at an extended blackbody source which is 100$ 
polarized , such that the temperature in the direction of polarization 
is equal to the temperature of the eccosorb and such that the 
temperature in the orthogonal direction is equal to the temperature 
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of the sky. If this source were aligned with the symmetry axis of the 

switch, it would generate a signal at the second detector at the 

switching frequency, like that in Figure 2.8. In this way one can say 

that a polarized source of magnitude T , -T , will produce a 

* eccosorb sky 

signal at the second detector whose effective first harmonic peak-to- 

peak value is I.l64 (V -V , ). where the calibration value 

e ' eccosorb sky 

derived in Appendix B2 has been used. Redefining C to accommodate 
the calibration factor 9 it can be written 




1.164 (v -v , ) 

v sky eccosorb 

(T -T 7) 

sky eccosorb 



( 2 . 22 ) 



where a typical value for V , -V , was 7*0 mv measured to an 

sky eccosorb 

accuracy of .1 mv and for T , -T . was -290. 0°K measured to 

sky eccosoro 

an accuracy of ,25°K. 

Next, I insert a calibration signal at the second detector. This 
is a sinewave of known peak-to-peak amplitude, which is derived from 
the audio oscillator through a special circuit. With the lock-in 
phase adjusted for maximum signal, the output voltage is read and a 
second parameter which gives the gain of the rest of the receiver, 
is computed: 



C 2 V 



lock-in 



(2.23) 



calibration signal 
The product of and then gives the gain of the entire polarimeter. 

In practice, the measurement was done utilizing the output of the 
entire data collection system rather than just that of the lock-in. 

In this way the overall gain of the entire system could be measured 
and individual intermediate steps, where human reading errors could 
multiply, were minimized. 



SIN (1.84 cose 



36 




o 



FIG. 2.8 

OUTPUT WAVEFORM FROM SWITCH 
FOR MAXIMUM SIGNAL 
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Although the temperature of the eccosorb was always known, that of 
the sky couldn't be with great certainty. For this reason, I picked a 
standard temperature of 5°K for the sky for purposes of calibration, 
even though, typically, the actual temperature will fluctuate from 
5-10°K. It was felt that the error introduced by this assumption 
was only on the order of 5/300 or about 2$, which, compared to the 
uncertainty due to receiver noise alone, is negligible. 

After a few calibrations it became clear that the overall gain 
stability of the system was good enough that a calibration on the 
order of once a week was enough to ensure good accuracy. Table 2.1 
is a list of calibrations taken during the data collection period. 

The units in which is expressed are counts per second at the 

output of the voltage-to-frequency converter per millidegree Kelvin 
at the polarimeter input. The output of this device is directly 
proportional to the lock-in output and was used to make integration 
of the signal and digitization of the result much easier. The entries 
are segregated into two groups because work on the receiver changed 
the overall gain midway through the data taking period. Figure 2.9 
plots the distributions of Cg for the two groups, which appear 
approximately normally distributed. It is interesting to note that 
if the data had been interpreted using the mean of each group of 
calibrations, the expected error in the calibration, equal to the 
standard deviation of the measurements divided by the mean, would be 
about 2.8$. 

As a final comment, I would like to discuss the overall effective- 
ness of this method considering that there was no true polarized 
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TABLE 2.1 

POLARIMETER CALIBRATIONS 



First Group 








Date 


C C 2 (CNTS/mdeg K) 


Date ■ 


C 1 C 2 (CNTS/mdeg K) 


8/11/72 


40.50 


11/20/72 


40.31 


8/24/72 


39.91 


11/29/72 


40.91 


9/7/72 


42.28 


12/7/.72 


39.81 


9/18/72 


40.75 


12/14/72 


39-84 


9/26/72 


40.57 


12/29/72 


39-08 


10/2/72 


41.79 


1/4/73 


40.43 


10/10/72 


42.85 


1/7/73 


39.11 


10/23/72 


41.96 


1/11/73 


41.64 


10/30/72 


39.82 


1/18/73 


42.61 


11/9/72 


41.55 


1/26/73 


41.75 


C^C 2 = 40.83 


CT C 1 C 2 = 1.14 




2 = .028 



Second Group 



Date 



2/4/73 

2/13/73 

2/20/73 

2/28/73 

3/9/73 

3/20/73 

3/28/73 

4/6/73 

4/16/73 



(CNTS/mdeg K) 


Date 


C C 2 (CNTS/mdeg K) 


43.93 


4/24/73 


45.00 


44.42 


5/6/73 


44.49 


42.40 


5/29/73 


43.32 


42.72 


6/11/73 


44.81 


44.12 


7/4/73 


45.63 


44.21 


7/20/73 


44.07 


44.16 


8/1/73 


45.24 


44.48 


8/10/73 


47.47 


43.00 






ct c 1 c 2 = 1.19 


CT c n cV c“c 
1 2 1 


2 = 0,27 



C 1 C 2 = 44,32 
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DISTRIBUTION OF CALIBRATION 
CONSTANTS C< C 2 
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FIG. 2.9 
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standard against which to calibrate. As a means to justify the method, 
I made a comparison between the results given by my polarimeter and the 
results given by a polarization sensitive power radiometer, both 
measuring the same polarized source. Used as a source was a klystron 
radiating through a rectangular standard gain horn suspended four 

stories above the measuring devices and giving an almost fully 

« 

polarized signal. A reading with the polarimeter was sandwiched in 
between two readings taken with the power radiometer. Agreement 
between the two radiometer readings was 2 while agreement between 
their mean and the polarimeter reading was 4$. Because the power 
radiometer and the polarimeter used most of the same components (the 
only difference being the addition of a second dual mode transducer 
and comparison load between the antenna and the mouth of the switch), 
changing configuration meant partial dismantling of the apparatus. 
During this time great care was taken to maintain alignment with the 
source, but not all the microwave plumbing was uniform. For this 
reason it was felt that the errors involved were ones of alignment 
and not ones of calibration. The error in each measurement due to 
receiver noise was only about 15 mdeg.K, while the signal was about 
280 °K. Given this, the discrepancy between the two power radiometer 
measurements must also have been due to alignment, as the calibration 
constant only varied about .1$ between the two runs. 
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III . PROCEDURE 



3.1 EXPERIMENTAL PROCEDURE 

Eor the reasons explained in previous sections, data were taken 
by the polarimeter with its axis mounted vertically and observing the 
zenith. Rotation of the earth sweeps out a circle of constant 
declination in the sky. By rotating the poldrimeter by 45 degrees 
between observations, measurement of the Q and U Stokes parameters 
can be made. 

There is one additional source of systematic error inherent in 
this method. This source can be eliminated by data talcing procedures. 

In Chapter 2 and Appendix B2, I discuss at length the problems of 
external parameters influencing the offset of the Faraday rotation 
switch. Aside from these problems, it is easy to show that any 
offset at all is harmful. For example, for a switch offset of 50 
millidegrees entering the receiver, any gain fluctuations on the 
order of 2°j 0 will produce systematic effects of around a millidegree 
in the data, an unacceptable level. To remove this we must find some 
way of eliminating the offset from the data in such a way that the 
effects of nonrandom gain fluctuations will be eliminated. 

A possible solution utilizes the symmetry properties of the 
switch itself. With the polarimeter aligned with one polarization 
axis north-south, it produces a signal proportional to 

T kto -T„ tt +T „„ , . Rotating the axis of the polarimeter by ninety degrees 

jWo iiiW oiiset 

produces a signal T^- T^ s + where the contribution due to 

^'offset < ^ oesn '^ change. It is a systematic effect in the polarimeter 
itself, and is independent of external orientation if the switch is 
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properly shielded. The readings taken in the two positions are 
subtracted to give 




from which the offset has been eliminated. Of course, how well this 
subtraction works depends on how constant the offset is over the 

rotation timescale. Keeping the integration , times short reduces 

■> 

sensitivity to longer term drifts in the offset. 

Combining the idea of rotating the polarimeter by ninety degrees 
between equivalent readings with the need to take readings separated 
by forty-five degrees, one arrives at a natural solution to both 
problems; that is, by rotating by forty-five degree steps around the 
compass. In this way one measures both Stokes parameters and for 
each measurement there are two others separated from it by ninety degrees, 
to subtract out switch offset. 

The polarimeter was mounted on a rotating platform with magnetically 
actuated reed switches located at forty- five degree intervals around 
the periphery, and labeled one through eight. Data taking was con- 
ducted by rotating the table in a clockwise direction until it actuated 
a switch, stopping and integrating for five minutes at that position, 
and then recording the result and moving on to the next position. 

When the recording at position eight was completed, the table would 
reverse direction and return to position one. 

To obtain the most useful results, it was felt that the following 
experimental quantities had to be measured: 
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1. The output of the integrator, which sums the polarimeter 

output for five minutes e 

2. The time of the observation. This quantity can be easily 

converted to local mean sidereal time which is 
identical to the right ascension of the observation. 

3. The angular position of the polarimeter, which determines 

which Stokes parameter one is measuring and with which 
sign. 

4. The temperatures of the switch, the inside of the observa- 

tion shed, and the outside environment. The first 
of these was needed to correlate the switch offset 
with temperature. I also kept track of the inside 
shed temperature to insure that my temperature 
controller was working properly. In the event that 
ground radiation became a large problem, measuring 
the temperature of the surrounding buildings would 
help in removing the effect. 

The apparatus involved in taking and recording these observations 
is the subject of the next few sections. Figure 3.1 is a diagram of 
the entire system. It was felt that by automating as much as 
possible, the chances for human error could be reduced, especially 
in an undertaking of this type where the observations extend over 
many months. 



POLARIMETER 
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DATA COLLECTION SYSTEM 

FiG.3.1 
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3.2 PREPARATION OF SIGNAL 

The output of the lock-in amplifier was recorded in two ways. 
First the output itself was recorded directly so that any abnormal- 
ities in function or signal received could be detected and evaluated. 
Two examples which were easiest to detect in this way were rain and 
local thunderstorm activity. The most dependable way found for 
recording was a hand wound, galvanometer movement, Esterline Angus 
chart recorder which could be left unattended for two to three days 
at a time. Combined with a marking device which delineated the 
integration periods in each position and the hours of the day, I was 
able to relate a disturbance in the apparatus to a recorded data 
point . 

To integrate the output of the lock-in, I used a prescaler and 
voltage -to -frequency converter which changed a -10 to +10 volt 
output to a 0 to 100,000 Hz signal. Integration was performed by 
a remotely controlled scaler, which was capable of integrating the 
signal for periods from zero to 9999 seconds, could divide the input 
by powers of ten, and could display six digits of output information 
in binary coded decimal form. This result was presented in parallel 
form to a parallel to serial converter which prepared it for 



recording . 
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3.3 TIMING 

Time information was provided by a clock built by E. Groth and 
used primarily for pulsar timing. Needless to say its micro-second 
accuracy was not needed, but it did provide time in days, hours, 
minutes and seconds of Eastern Standard Time relative to epoch 
1970.0. Upon completion of each integration period, the output 
of the clock was stored and the result presented to the parallel - 
to-serial converter for recording. This provided a convenient 
and reproducable point from which to calculate the time of the 
middle of each integration run. 

3.4 TEMPERA rURE MEASUREMENT 

The temperature at key locations was measured using four 
linearized thermistor composites connected to a measuring circuit 
and digital voltmeter through a multiplexing scanner. During each 
run one of the thermistor probes was measured and recorded. The 
output of the digital voltmeter and the position of the scanner were 
presented to the parallel-to-serial converter for recording. It was 
thought that this information might be used to correlate out the 
temperature dependent effects . 
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3.5 RECORDING OF DATA 

The data outputs are ultimately recorded in IBM compatible form 
on digital tape. The intermediary in the process is the parallel-to- 
serial converter. When strobed, it takes all the data presented to it 
in parallel form and feeds it serially one byte at a time to an 
incremental nine track 800 tyte per inch tape unit. In addition, it 
controls the blocking of the data into files knd records to facilitate 
analysis. This device plus the scaler and dvm scanner were constructed 
by Karl C. Davis. 

3.6 EXPERIMENT CONTROL 

Overall control and sequencing of the experiment was performed by 
an experiment control box, the function of which was to command the 
polarimeter, the integrator and parallel-to-serial converter and 
sequence their functions properly. The order of operations while the 
experiment was being conducted was : 

1. Rotate to a new position. If at position eight, rewind 

until reaching position eight again and then step 
forward to position one. 

2. Start integration. 

3. Stop integration and record data. 

4. Recording finished, rotate to a new position. 

Cycle time for the apparatus with a three hundred second integration 
time was three hundred and three seconds except when rewinding, in which 
case it was three hundred and twenty seconds. The device would maintain 
this timing without losing a second for a week or more, making it very 
easy to spot malfunctions and timing errors. 
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All switching in the system is either low voltage or solid state 
to avoid radio frequency interference with the sensitive receiver. 

For this reason, AC switching for the polarimeter rotation motors is 
tied to the zero crossing of the line voltage, using a semi- 






conductor triac. 
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IV. ANALYSIS AND INTERPRETATION 



4.1 DATA EDITING 

Data was taken for this experiment from 1 February 1972 until 
10 August 1973. The data before 9 August 1972 was mainly for test and 
evaluation of the device and has been discarded out of hand. The data 
since that time is considered suitable for analysis and has been 
searched for the effects of cosmological polarization. The twelve 
month data set contains approximately 94,000 individual polarization 
measurements, of which about half were devoted to each Stokes 
parameter. As in any experiment of this length, bad weather, equip- 
ment malfunctions, and human errors took their toll, resulting in 
discarded readings. 

My guiding principle in throwing away data was not to, if at all 
possible. When necessary it would have been best to completely auto- 
mate the procedure in order to avoid human intervention with its 
attendant bias. But at certain points, in particular the screening 
of chart recordings of the lock-in output, there was.no other way. 

In general there were four sources of abnormality in the lock-in 
output : 

1. The most common and troublesome was rain. The polarimeter 

was insensitive to clouds and other forms of bad weather, 
but precipitation would collect on the mylar window covering 
the apparatus and appear as a strongly polarized signal. 

This would many times drive the lock-in off scale. 

Data showing this type of behavior was rejected along 
with enough data on either side to insure that all 



50 



contamination was removed. 

2. The next most frequent source of interference was thunder- 

storm activity, especially in the summer months. This 
was characterized by strong and repeated noise spikes 
in the output signal, which caused' the data to be 
discarded. Isolated noise spikes were not thrown 
away by manual means. Later stages of automated 
editing evaluated their effect and provided more 
impartial criteria for their elimination. 

3. During the summer months when the sun reaches altitudes 

close to the zenith, its presence as a hot source 
produces a lock-in signal when acting on the non- 
symmetric antenna gain pattern. Again this data 
was not deleted manually but was handled later when 
all data with the sun above the ground shield and 
geometrically illuminating the antenna were discarded 

4. Finally, there would be an occasional malfunction in the 

receiver and its associated equipment, which would 
usually be indicated by a lack of receiver noise in 
the output trace. This was easily noticed, the 
receiver repaired and the data removed. 

A miscellaneous problem that does not fit into any of the above 
categories, but which at one point had me convinced, that the apparatus 
was broken, was that of a small very inconspicuous spider. By slowly 
crawling around in the antenna throat he did a great imitation of 
serious receiver gain drifts. Naturally his virtuoso performance 



was eliminated. 
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Removal of bad data constituted the first stage of editing. Next 
I searched for timing errors to make sure that the automated data 
collection apparatus was functioning properly and that the integration 
periods were uniform. These events fell into two categories, clock 
errors and actual erratic operation. In the former case the equipment 
would be operating normally, but the clock would indicate the wrong 
time. Resetting of the clock for pulsar timing was the most common 
cause of this and would appear as a step function in the cycle, 
which would be followed by a step with the opposite sign when the 
correction was completed. The cycle time of the apparatus was 
consistent enough to enable one to reconstruct the correct time 
even when the step remained for several hours. 

Harder to correct was difficulty due to erratic behavior of the 
apparatus, which was usually associated with thunderstorm activity 
or other forms of electrical interference. In most of these cases, 
where the orderly sequence of readings could not be restored, the 
errors were left for a later stage where they would automatically 
be eliminated. 

On the average of every one or two weeks, the cycle would lose a 
second or two in phase for no apparent reason. No attempt was made 
to correct for this. In later processing, isolated errors in timing 
of up to five seconds were allowed to pass. A change in the 
integration period of up to five seconds means an error in the 
measurement of under two percent or far less than the uncertainty 



due to receiver noise alone. 
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As a final editing procedure, the problem of isolated interference 
spikes in the data was considered. Preliminary examination showed that 
after removal of switch offset, the averaged signal would be within 
10 millidegrees of zero. The minimum detectable signal or standard 
deviation of the noise for a single processed reading was about 15 
millidegrees based on measurements of scatter in the data. The 
clipping level for rejection of spikes was set three full standard 
deviations above the 10 millidegree level, at 56 millidegrees Kelvin. 
Only 387 points had to be eliminated in this way. 

In all, out of about 9^000 points collected, only 38,029 points 
or about kO .4$ of the data were actually used. Unfortunately they 
were not uniformly distributed throughout the year for two reasons. 
First, more data get thrown out in summer when the sun is higher in 
the sky. Worsening the effect was a very wet spring and early summer 
in the Princeton area, combined with unstable air and much electrical 



storm activity. 
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4.2 PRELIMINARY PROCESSING 

The final stage of data preparation is carrying out the subtraction 
(described under experimental procedure) for removal of the switch 
offset. At this point we have two roughly equal sets of measurements 
of the two Stokes parameters describing the polarization. Each of 
these sets contains measurements at four of the polarimeter positions. 
One set comprises all the odd numbered switch positions and the other 
has the even, as in Figure 4.1. For the Q parameters, that position 
closest to North or position two, was called positive in accordance 
with normal convention. It was found by checking with a polarized 
source placed due North of the apparatus that this position actually 
pointed 11.75 degrees to the Northeast of North. I will still call 
this the Q, ' parameter (primed to distinguish it) and correct the 
discrepancy at a later stage. For the U' parameters, the positive 
direction is that closest to the Northeast or position three. 

Positions 180 degrees relative to two and three have the same sign 
by symmetry, while those 90 degrees away have the opposite sign. 

To perform the subtraction, each reading was assigned its 
characteristic sign and then adjacent pairs of measurements within 
each set were averaged as in Figure 4.2. The time of observation 
assigned to each averaged pair is the average c£ the midpoints of 
each integration period. 

If we look at this operation from another point of view, we see 
that what we are doing is multiplying the data by a periodic signal 
of frequency twice that of the polarimeter rotation frequency. This 
is exactly what is done in the mixer stage of a lock-in amplifier. 



NORTH 




POLARIMETER POSITIONS 
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OLARIMETER 
POSITION 



1 

2 

3 

4 

5 

6 

7 

8 
1 
2 

3 

4 

5 

6 

7 

8 



INTEGRATOR 



TIME 

15=45=39 


OUTPUT 

1628.97 


Q ' 


U' 


15=50=42 


1663.33 




-1.54 


15=55=45 


1625.90 


0.33 




16=00=48 


1663.98 




-23.67 


16=55=51 


1673.23 


1.13 




16=10=54 


1661.72 




11.31 


16=15 = 57 


1695.85 


7.89 




16=21=00 


1677.49 




-8.43 


16=26=20 


1712.71 


12.33 




16=31=23 


1652.83 




-25.89 


16=36=26 


1660.93 


-6.06 




16=41 = 29 


1640.70 




4.84 


16=46=32 


1651.24 


-21.76 




16=51 = 35 


1684.21 




1.69 


16=56=38 


1654.63 


26.57 




17=01=41 


1737.35 







TIME 

OF 

AVERAGE 

15=48 = 12 
15 = 53=15 

15 = 58=18 
16=03=21 
16=08=24 
16= 13=27 

16=18=38 
16= 23 = 41 
16=28 = 53 
16=33=56 
16=38 = 59 

16 = 44 = 02 
16=49 = 05 
16=54=08 



SUBTRACTION SCHEME 

FIG. 4.2 
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Because we do the Q,' and U' parameters separately, we extract both 
phases of the signal simultaneously just as in a two phase or vector 
lock-in. To visualize things a little more clearly, consider the 
effect on incident DC polarized radiation as the polarimeter is 
rotated. The rotation causes a sinusoidal variation in intensity. 
Because of the 180 degree symmetry of the switch, this signal has twice 
the rotation frequency of the polarimeter. JIow by multiplying by a 
periodic signal the way we have done, we demodulate the signal and 
extract the fundamental Fourier component for each phase, exactly as 
is done in a lock-in. A double lock-in technique of this type gives 
good immunity against effects from gain fluctuations and drifts in 
the switch offset } which are on a time scale longer than half the 
polarimeter rotation period, or twenty minutes. 

4.3 ANALYSIS 

I now have two sets of independent measurements, the Q' and U' 
parameters, each with a time resolution of about ten minutes in 
sidereal time. One might think that this could be interpreted as a 
spatial resolution in terms of right ascension. This is not true 
because of the low resolution of the antenna, which is only about one 
hour in right ascension. To bin the data on any sidereal time scale 
shorter than this is to attempt to display resolution which does not 
exist. The next logical step is to bin the data at its true resolution 
limit, one hour in right ascension (local sidereal time). 

In ord.er to check for serious solar effects, the binning was also 
carried out in solar time. After binning, Q,' and U' were transformed 
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to the true Q and U parameters before further analysis. The results 
of the binning are shown in Figure 4.3 for solar time and Figure 4.4 
for sidereal time. The temporal or spatial center of each bin may be 
obtained by subtracting thirty minutes from the bin number. In other 
words , bin four contains all the data between three and four hours . 

In each set of data there is a DC level, which when combined from 
both sets gives a resultant DC polarization yector with magnitude of 
1.10 ± .14 millidegrees Kelvin and direction 116.34 ± 3.64 degrees from 
North. On the reciprocal of this bearing is the most prominent feature 
on the horizon, a tower containing elevator motors on the roof of Jadwin 
Physics Laboratory. VJhen reviewing the data in smaller segments, it was 
found that this level wasn't constant throughout the year of observation. 
During the first two months, the DC vector was 6.6 ± .28 millidegrees with 
a direction of 150.46 ±1.21 degrees. Other than this, no other abnormal- 
ities were found. Closer examination showed that the change took place 
sometime close to the tenth of October. Checking my notebook, I dis- 
covered that on the fifteenth of October I had decided that the calibra- 
tion probe installed in the antenna wasn't being used and should be removed. 
It appeared that the presence of a stub in the horn had perturbed the 
gain pattern enough to significantly alter the response to ground 
radiation, in effect, making the polarimeter six times more sensitive 
to the elevator tower radiation. The discrepancy in angle of about 
+36 degrees can be traced to the fact that the calibration probe was 
located approximately -40 degrees from the positive symmetry plane of 

the switch. VJhen the polarimeter was aligned at 150 degrees, the 

/ 

probe was aimed at the tower. A check was performed with a slotted 



m °K 
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3.2 POLARIZATION DATA 
FOLDED INTO SOLAR BINS 
BIN WIDTH = 1 hr 




BIN 



FIG. 4.3 
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3.2 POLARIZATION DATA 
FOLDED INTO SIDEREAL BINS 



BIN WIDTH = 1 hr 




BIN 

FIG. 4.4 
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plate providing a polarized source to determine whether or not the probe 
changed the effective position angle of the switch's symmetry plane. 

No change could be detected between the position angle with or without 
the probe. 

Because the only effect of the probe was to change the DC level 
of the first two months data, this effect was subtracted out by 
adjusting the level of these data to agree wjLth the average level of 
the rest of the year. Failure to do this would have left a step 
function in the data and produced possible harmonic contamination from 
the spectrum of frequencies contained in a step. The only other 
alternative would have been to discard the data and no other reason 
to do this could be found. 

It seems fairly certain that a good portion of the DC signal is 
due to ground radiation. From equation 2.l6 with (3 equal to zero we 
obtain the expression for the Q Stokes parameter, given an axially 
symmetric universe. This parameter has a DC component given by 



where 6 is the declination and § is the angle between the rotation 
axis of the earth and the symmetry axis of the universe. For (3 equal 
to 45 degrees we get the U Stokes parameter for which the predicted 
DC value is zero. The experimentally measured values for these two 
quantities are 




(4.1) 



Q dc = - .67 ± .14 mdeg K, 
UL_ = - .88 dh .14 mdeg K. 

JA/ 



(4.2) 
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It's interesting to note that the U component which is supposed to be 
zero is larger than Q, the non-zero component predicted by the model. 
Taking this with the fact that the resulting vector points right at 
the elevator tower argues for this being a local effect due to ground 
radiation. Yet the possibility of a nonzero contribution from the 
galaxy or an asymmetrically expanding universe cannot be ruled out. 

The estimates of the ground radiation are nowhere near accurate enough 
to allow separating out the ground contribution in an effect this 
small. For this reason the result for Q, in(4.e) must be taken as an 
upper limit on the DC polarization from an asymmetrically expanding 
universe: 

Q dc < .85 mdeg K, with 90 $ confidence. 

This ends the analysis of the DC part of the signal. Next I 
will subtract the respective DC component from each Q and U value 
and investigate the result for evidence of a cosmological polarization. 

The result of the subtraction is given in Table 4.1. 

Now I must ask whether or not there is any significance to what 
is left after the DC contribution is subtracted? This poses a statistical 
question. To be answered, the statistics of the data must be investigated. 
The Q, and U parameters in Table 4.1 are averages of samples drawn from 
normally distributed populations whose means are the true signal in 
each bin. Because the samples are large (on the order of several 
hundred members), the population variance can be replaced by the sample 
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variance defined as 
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TABLE 4.1 

STOKES PARAMETERS AFTER SUBTRACTION OF DC POLARIZATION 



Q 



T - T 
NS EW 



o = Experimental error in Q, 
Q 



U = T -T 

RE,SW NW,SE 



Oy = Experimental error in U 



N~ = Number of measurements of Q, in average 

0 , 

Ny = Number of measurements of U in average 



Bin 

Folded 


Q °Q 

in Solar Bins 


n q 


U 


°u 


N u 


1 


O.98 


0.61 


1034 


-O.19 


0.60 


1035 


2 


-o.4o 


0.6l 


1044 


-0.31 


0.61 


1047 


3 


0.16 


0.60 


1051 


-0.24 


0.60 


1054 


4 


-0.24 


0.57 


1053 


0.62 


0.59 


1055 


5 


1.27 


0.58 


1050 


0.22 


0.58 


1050 


6 


0.16 


0.59 


1055 


-0.51 


0.59 


1048 


7 


-0.54 


O.63 


1010 


0.52 


0.62 


1012 


8 


-0.94 


O.65 


936 


0.54 


O.65 


928 


9 


1.32 


0.72 


635 


-0.71 


0.67 


636 


10 


-0.33 


0.87 


445 


0.45 


• O.85 


447 


ll 


-0.46 


1.10 


283 


-1.84 


1.10 


279 


12 


-0.57 


1.27 


224 


-1.07 


I.27 


217 


13 


-0.46 


1.24 


211 


-0.13 


1.24 


2l4 


14 


-0.53 


1.10 


274 


-0.73 


l.ll 


274 


15 


-1.82 


0.93 


437 


-O.63 


0.95 


436 


l6 


-0.68 


0.80 


578 


0.01 


0.79 


568 


17 


-0.54 


0.69 


830 


O.89 


0.10 


833 


18 


0.82 


0.71 


936 


0.10 


0.65 


933 


19 


-0.73 


0.67 


933 


0.25 


0.67 


918 


20 


0.10 


0.64 


963 


0.57 


0.65 


962 


21 


0.90 


0.61 


975 


-1.34 


0.61 


975 


22 


0.16 


0.60 


1000 


o.4o 


0.60 


1003 


23 


0.72 


0.60 


1037 


0.33 


o.59 


1034 


24 


-1.25 


0.62 


1038 


O.16 


0.60 


1039 
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TABLE 4.1 (continued) 



Bin 


Q 


°q 


N- 

Q, 


U 


°u 


N u 


Folded 


in Sidereal Bins 




1 


0.0 


0.67 


779 


O.67 


0.66 


776 


2 


- 0.24 


0.67 


813 


0.24 


0.65 


805 


3 


0.44 


0.72 


813 


1.53 


0.68 


811 


4 


-1.62 


0.66 


809 


1.13 


0.67 


813 


5 


0.65 


0.65 


810 


0.47 


0.67 


815 


6 


- 0.13 


0.64 


809 


0.98 


0.67 


800 


7 


- 0.15 


0.66 


769 


0.21 


0.66 


773 


8 


- 1.40 


0.70 


730 


0.06 


O.67 


718 


9 


1.14 


0.72 


706 


- 0.17 


0.74 


704 


10 


-0.09 


0.74 


681 


0.03 


0.74 


681 


11 


0.83 


0.74 


664 


-0.26 


0.72 


663 


12 


0.09 


0.73 


648 


-O.29 


0.74 


64 l 


13 


1.06 


0.80 


646 


0.04 


0.80 


64 l 


14 


0.75 


0.82 


643 


0.72 


0.80 


642 


15 


1.89 


0.78 


708 


- 0.64 


0.77 


710 


l6 


O.19 


0.68 


838 


- 0.95 


0.68 


847 


17 


0 .i 4 


0.62 


1003 


-0.56 


0.62 


989 


18 


-0.23 


0.62 


1024 


-0.08 


0.62 


1014 


19 


-0.87 


0.60 


1022 


0.18 


0.58 


1021 


20 


0.15 


0.63 


1008 


-O.69 


0.63 


1027 


21 


1.10 


0.69 


842 


- 0.01 


0.71 


839 


22 


- 0.55 


0.72 


762 


-1.66 


0.72 


773 


23 


- 0.79 


0.73 


732 


- 1.45 


0.73 


734 


24 


-I.03 


0.70 


773 


o.i 4 


0.70 


760 
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N 



— N 2 



= A - x > 

~ I 1 < x i - x > 2 



(4.4) 



It should be noted that this is not the variance of the original 
measurements and is therefore not a fair estimate of the experimental 
error. When I performed the subtraction between points 90 degrees 
apart in position in order to subtract out the switch offset, I formed 
a new data set by averaging adjacent points in the old. When each 
point is assigned its characteristic sign dependent upon polarimeter 
position, the process may be written as 

T ± = P i + F i + 1, (4.5) 

2 

where P is a given Stokes parameter. By the propagation of errors one 
can show that the variance of the new data set is related to the old 
by 

2 °P 2 

a t = JL __ • (4.6) 

i 2 

Thus in terms of the original measurements , Equation 4.4 is an estimate 
of the variance divided by two; therefore , I must multiply by two 
to get back a true estimate of the scatter in the data. 

There still remains one step left in order to determine the error 
in the measurements as they are binned. As the data is folded, five 
consecutive readings from each day's data are averaged into a given bin. 
Because of the method of subtraction, these five points are derived 
from six independent points in the original set of measurements. 

Because a whole day intervenes between each successive group of five 
measurements in a particular bin, these sets of five are independent 
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and averaging over all the points in a bin is just equivalent to averaging 
together the averages of each day's contribution. This daily average 
can be written as 



P.+ 2 P. .+ 2 P. 0 + 2 P. 0 + 2 P. |,+ P. _ 

, _ _i l+l i+2 i+3 1+4 i+5 

F Avg " 10 ■ ’ 



( 4 . 7 ) 



and has a variance given by 



a P' 



P. 

l 



(4.8) 



Avg 556 



where Op is the variance of the original measurements. Just computing 
i 

Equation 4.4 for the bin in question and dividing by the number of 



points would give p 



P! 

Avg 



( 4 . 9 ) 



Therefore, in order to compute the variance of the mean for each of my 

bins, I must compute Op, from Equation 4.4, multiply by 2/N to get 
2 1 

o_ , and finally multiply by 5 . 0 / 5*56 to correct for the binning 
P. 

procedure . 

In fitting a DC level to each set of Q, and U, I used the method of 
least squares, which is to say, I minimized the function 



24 (P. - P) 2 
F(P) = E ^73 

i=1 CT ?. 

1 



(4.10) 



where o_ refers to the corrected variance of the mean described above. 
P. 



1 

If P is the true mean of the population in each bin, then F(P) is 
2 

distributed as x with 24-1-23 degrees of freedom. The hypothesis 
that P is the mean of each bin is just the statement that the polariza- 
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tion is isotropic for that particular Stokes parameter. Comparing the 
actual minimized value of F(p) to the x distribution function for 
23 degrees, tests whether the isotropic hypothesis is valid. The 
values of F(P)^ for each of the fits to the Q and U Stokes para- 
meters and the probabilities that each value or greater would be 
given by a set of random numbers are : 

SOLAR BINS SIDEREAL BINS 





2 

X 


Probability 


2 

X 


Probability 


Q: 


30.60 




32.98 


10 % 


U: 


17.93 


1 % 


27.07 


2 % 



Except for the U parameter binned in solar time, the probabilities of 
the parameters having a completely random character is small. This 
means that statistically the measured polarization is not isotropic, 
or at least the probability of its being isotropic on the basis of 
this data is small. I must now determine whether this anisotropy 
has any real meaning for cosmology or whether it can be explained in 
terms of more local effects. In our model of an axially symmetric 
universe, in addition to the DC polarization, we are interested in the 
sinusoidal variations of the Q, and U parameters which have 12 and 24 
hour sidereal periods. Because I am dealing with a non-uniform, finite 
data string, I must also be concerned with contamination of the 
sidereal contributions by strong signals in the corresponding solar 
frequencies . 

The last effect can arise in two ways. First, because we are 
dealing with finite data sets, adjacent bins in the power spectrum 
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will not be completely uncorrelated. A strong peak in one bin will 
enhance an adjacent bin by as much as 15 percent of the value of the 
strong peak. This occurs because the elementary frequency bins are 
no longer rectangular in shape, but fall off as sin x/x where x is 
in units of one over the length of the data string in seconds . When 
the data sets are not uniform in variance over their length (as when 
data are thrown away), the sine and cosine functions used to extract 
the Fourier components from which the power spectrum is derived are no 
longer orthogonal. This increases the correlation between adjacent 
bins. For this reason we must be careful that all large solar effects 
are screened out. This is particularly true of the 24 hour solar 
and sidereal components, which are adjacent elementary bins in the 
power spectrum after one year of data. 

Another way in which solar effects can hurt is through seasonal 
effects modulating the daily solar cycles. The product of a 24 month 
seasonal variation with a solar period gives a beat frequency equal 
to the corresponding sidereal period. To look for these effects, 
the best starting point is to examine the data folded into solar 
bins for a signal which would be large enough to contaminate the 
sidereal results. 

Table 4.2 gives the results of fitting the eight lowest 
frequency sine waves to the data in both solar and sidereal time. 

The procedure was to minimize the function 

p 

24 (P^- A^Sinootw - B^cosut^) 
ill ~~2 ’ 

a P. 

1 
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with respect to the parameters A^ and B^. From A^ and one can 
compute an amplitude and a phase for the best fit sine function 






Amplitude = (A f + B f ) 

Phase = ARCTAN (^f ) , 
A f 



(4.12) 



for each frequency component. Starting with the 24 hour component, each 

wave was computed, subtracted from the data, and the residuals used to 
2 2 

compute ax. At each step I compared the x' to a random distribution 
with two less degrees of freedom. The next fit then uses the residuals 
of the preceding step. Table 4.2 also gives the expected value for each 
frequency component. Because of the non-zero variance of each data point, 
there is an expected value of the amplitude of each frequency component 
given by 



<Amplitude>„ =|ig 



2 o 

4 sin out,. r ^4 cos wt. ^2 _ ^ £4 coswt.. sinou 



_i [(. . 
2 1=1 



i=J 



2 



in 



^4 cosout. 2,4 sinuut. 



'i=l 



a. 2 
x 



i=l 



~2~ 



,2,4 cosujt. sincut. x 2 n 2 

- 1=1 1= i > 3 



a. 2 

l 



2,4 cos cut. 2,4 sincut. 0 

& 



a. ‘ 

l 



a. ' 
i 



24 cosout . sincut. 0 

1 ) 3 



r 24 cosout. 24 sincut. 
L£ 1 £ _ i 

i=l p — i=l 



a. c 
i 



a. 

i 



2” 



,2,4 coscut . sincut. x 2n2 

<?=! i, i > J 

a. ^ 



(4.13) 



To test the statistical significance of an amplitude , one takes the 
difference between the amplitude and its expected value for random data 
and compares this with the error in the amplitude. 

Investigation of the solar amplitudes shows only one with a 
significant value, the 4 hour period in the Q parameter. This 
frequency is well separated from its sidereal equivalent in the power 

spectrum, as one needs only two months of data to resolve the two. 
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TABLE 4.2 

HARMONIC ANALYSIS OP DATA 



P 

T 

CT T 



v 

2 

X 

<T> 



Period in hours of solar or sidereal time. 

Amplitude of harmonic component in mdeg K. 

Error in amplitude. 

Phase of component in degrees from Bin 0 (Bin 24). 

Error in phase. 

Degrees of freedom in fit. 

Chi-square for the fit. 

Expected value of T based on experimental errors and random 
numbers . 



T-<T> 

= Number of sigma of amplitude above expected value 

J- 



p 


T 


°T 


§ 


cr, 

$ 


V 


X 


<T> 


T-<T> 
















°T 


Q in 


Solar 


Bins 












24 


0.33 


0.21 


97.36 


34.35 


21 


26.94 


0.29 


0.19 


12 


0.19 


0.20 


-123.61 


62.45 


19 


25.59 


0.29 


-0.50 


8 


0.28 


0.20 


-42.31 


42.03 


17 


24.21 


0.29 


-0.05 


6 


0.22 


0.20 


-100 . 89 


53.08 


15 


23.25 


0.29 


-0.35 


4.8 


0.09 


0.20 


182.42 


129.61 


13 


23.22 ‘ 


0.29 


-1.00 


4 


0.60 


0.20 


73.14 


19.40 


11 


14.95 


0.29 


1.55 


3.4 


0.23 


0.20 


30.88 


51.09 


9 


14.96 


0.29 


-0.30 


3 


0.15 


0.20 


-89.83 


78.45 


7 


14.64 


0.29 


-0.70 


U in 


Solar 


Bins 














24 


0.10 


0.21 


123.20 


109.79 


21 


16.99 


0.29 


-0.90 


12 


0.23 


0.20 


-40.31 


49.17 


19 


15.35 


0.29 


-0.30 


8 


0.14 


0.20 


139-9^ 


84.06 


17 


14.45 


0.29 


-0.75 


6 


0.19 


0.20 


-163.76 


60.61 


15 


13.92 


0.29 


-0.50 


4.8 


0.31 


0.20 


-102.86 


37.02 


13 


11.79 


0.29 


0.10 


4 


0.24 


0.20 


-133.50 


48.93 


11 


10.96 


0.29 


-0.25 


3.4 


0.08 


0.20 


-124.28 


143.00 


9 


IO.92 


0.29 


-1.05 


3 


0.47 


0.20 


52.36 


24.81 


7 


5.61 


0.29 


0.90 
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TABLE 4.2 ( C ont inued ) 



p 


T 


°T 


§ 




V 


2 

X 


<T> 


T-<T> 


Q in 


Sidereal 


Bins 












°T 


24 


0.52 


0.21 


-80.52 


21.79 


21 


25.65 


0.29 


1.10 


12 


0.20 


0.20 


89.61 


58.59 


19 


24.35 


0.29 


-0.45 


8 


0.19 


0.20 


-83.75 


61.66 


17 


23.10 


0.29 


-0.50 


6 


0.31 


0.20 


-16.64 


37.81 


15 


20.89 


0.29 


0.10 


4.8 


0.4l 


0.20 


84.80 


28.58 


13 


17.11 


0.29 


0.60 


4 


0.36 


0.20 


62.64 


31.95 


11 


13.68 


0.29 


0.35 


3.4 


0.10 


0.20 


-59.81 


116.81 


9 


13.63 


0.29 


-0.95 


3 


0.39 


0.20 


-105.69 


30.13 


7 


10.42 


0.29 


0.50 


U in 


Sidereal 


Bins 














24 


0.58 


0.20 


39-24 


20.53 


21 


17.30 


0.29 


1.45 


12 


0.45 


0.21 


22.62 


25.30 


19 


12.78 


0.29 


0.76 


8 


0.33 


0.20 


-1.32 


35.01 


17 


9.63 


0.29 


0.20 


6 


0.30 


0.20 


79.77 


38.67 


15 


7.14 


0.29 


0.05 


4.8 


0.17 


0.20 


187.52 


67.42 


13 


6.46 


0.29 


-0.60 


4 


0.08 


0.20 


40.47 


136.89 


11 


6.28 


0.29 


-1.05 


3.4 


0.29 


0.20 


183.51 


40.62 


9 


3.85 


0.29 


0.00 


3 


0.25 


0.20 


-163.63 


46.74 


7 


2.61 


0.29 


-0.20 
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The absence of significant solar values with frequencies lower than this 

indicates that contamination of sidereal frequencies with solar signal 

won T t be a problem. The solar signal is not zero as evidenced by the 
2 

high initial x but I would be a little surprised if it were, since 
most environmental variations which can affect the polarimeter are 
solar in nature. It appears that by shielding the polarimeter from 
ground radiation and by throwing away data when the sun was above the 
ground shield, I was able to avoid a large amount of contamination. 

Before discussing the data further, there is one feature of a fit 
of this type that should be noted. Ordinarily when we fit sines and 
cosines to a set of numbers with uniform variance, the sines and cosines 
are orthogonal functions over the set. In my case the variances of 
the points are not equal but depend on the number of measurements per 
point and the actual scatter in those measurements. To insure that I 
wasn't adding any systematics to the data in the fitting process, I 
compared the results obtained above to the results one would obtain 
assuming a constant variance across the data set. For each set I chose 
as a variance the average variance per point of the . corresponding 
Stokes parameter folded into sidereal bins. The results are given in 
Table 4 . 3 . The most notable differences are in the solar frequencies. 
Here the points closest to noon, which were deweighted because of 
fewer measurements (due to solar contamination), were given equal 
weight. It is also these points which are most affected by the sun 
and therefore, it is not surprising that the solar amplitudes become 
somewhat larger. The effect on the sidereal data is small. In both 
cases the Chi-squares seem to follow the same general progression as 
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TABLE 4.3 

HARMONIC ANALYSIS OF DATA ASSUMING UNIFORM VARIANCE 
P T $ v 

Q in Solar Bins 



24 


0.43 


93.32 


21 


27.44 


12 


0.31 


-114.01 


19 


24.97 


8 


0.22 


-13.12 


17 


23.68 


6 


0.17 


-109.24 


15 


22.92 


4.8 


0.12 


-117.86 t 


13 


22.58 


4 


0.59 


75.56 


11 


13.73 


3.4 


0.i4 


9.71 


9 


13.20 


3 


0.12 


-93.50 


7 


12.86 



U in Solar Bins 



24 


0.29 


121.77 


21 


21.17 


12 


0.39 


-46.24 


19 


17.26 


8 


0.24 


i4o.o8 


17 


15.72 


6 


0.08 


-164.68 


15 


15.55 


4.8 


0.33 


-105.00 


13 


12.71 


4 


0.12 


181.05 


11 


12.31 


3.4 


0.21 


-93.80 


9 


11.15 


3 


0.45 


66.11 


7 


5.96 


in Sidereal Bins 


24 


0.56 


-83.84 


21 


27.38 


12 


0.25 


82.76 


19 


25.78 


8 


0.22 


-83.32 


17 


24.58 


6 


0.34 


-8.13 


15 


21.66 


4.8 


0.39 


92.11 


13 


17.86 


4 


0.35 


56.00 


ll 


14.78 


3.4 


0.08 


-50.59 


9 


14.63 


3 


0.38 


-111.56 


7 


10.85 


in Sidereal Bins 


24 


0.59 


34.81 


21 


19.29 


12 


0.46 


20.38 


19 


13.68 


8 


0.35 


0.05 


17 


10.55 


6 


0.34 


78.64 


15 


7.63 


4.8 


0.18 


184.40 


13 


6.8l 


4 


0.10 


61.80 


ll 


6.56 


3.4 


0.34 


189.77 


9 


3.63 


3 


0.21 


-157-64 


7 


2.46 
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in the weighted analysis, indicating that no systematics were introduced 
due to the non-uniform weighting. 

Turning to the sidereal frequencies, there are three significant 
amplitudes, both 24 hour components and the 12 hour component for the 
U parameter. The largest of these is the 24 hour component of U. The 
positive maximum of the fitted sine wave falls at a right ascension of 
2.88 ± 1.37 hours which is about two hours a*?ay from the plane of the 
galaxy in the anticenter region. The negative minimum falls in the 
region of the Worth Galactic spur at 14.88 ± 1.37 hours. Looking at 
the surveys of galactic continuum polarization at lower frequencies, 
these regions seem to exhibit significant amounts of polarization. 
Looking also at the plot of the U data in Figure 4.4, we note that in 
addition to the high values around bins 3 and 4 and the low values 
around bins 15 and l6, there is another low region extending from bin 
20 to bin 23, or 19-5 to 22.5 hours in right ascension. Again looking 
at a galactic map, I find that the Cygnus arm of the galaxy, the 
hottest region to pass through my beam, falls at 20.5 hours. Because 
I cut the plane of the galaxy in two places roughly. 12 hours apart 
in right ascension, it isn't too difficult to imagine a 12 or 24 hour 
component being generated if any polarized component existed in these 
regions . 

49 

From a galactic polarization survey at 21 cm I found that 
individual measurements of 100 to 200 millidegrees were obtained in the 
regions I cover. Applying a spectral index of a -2.7 to this intensity, 
in order to extrapolate the reading to 3.2 cm, results in a possible 
polarization temperature of from .6 to 1.2 millidegrees Kelvin. 
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Although this arm waving proves little, it did. convince me that I 
had better look into polarized galactic radiation in detail. 



4.4 SUBTRACTION OF THE GALAXY CONTRIBUTION 

In any type of extrapolation, one would like to use data which has 
parameters as close to those of the data being corrected as possible. 
For a radio survey of the galaxy, the two most important parameters 
are frequency and angular resolution. Of these, frequency has the 
edge in importance when discussing a polarization survey. 

There are two ways in which frequency enters into extrapolating 
polarization measurements. Most directly it must be considered in 
terms of the spectrum of the radiation one is measuring. Most of the 
galactic polarized emission is thought to be non-thermal synchrotron 
radiation with a spectral index of -2.5 to -2.9. ? By spectral 

index we mean the quantity |3 in the expression 



T v 

£ - (r 1 ) P 



(4.14) 



1 1 

for scaling the temperature of radiation at different frequencies. If 
the radiation follows this simple relation, knowing the spectral index 
determines the magnitude of the polarization at one frequency, given 
the magnitude at another. 

The other aspect of polarized emission in the galaxy, the angle, 
also depends on frequency through the mechanism of Faraday rotation. 
The electron concentration in the galaxy together with the galactic 



magnetic field act to produce rotation of linear polarization in much 
the same way a magnetic field acting on a ferrite produces rotation 
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in the Faraday rotation switch. This effect is frequency dependent, 

53 

satisfying a relation of the form 

* = t + R a 2 , (4.15) 

Y Y o m 5 

where R is called the rotation measure and is given in units of 
m 

radians per meter squared, i|r and i|f are the angles of the radiation 

before and after rotation. In terms of the ^number of electrons and the 

longitudinal magnetic field along the line of sight, the rotation 

53 

measure may be written as 

R m = 8.1 x 10 5 J N B L dL RADIAiIS/m 2 . (4.l6) 

Because N and B are not constant with direction but vary in 
L 

a complex way, one cannot use them as fitting parameters to align two 

surveys at different frequencies. One needs at least two surveys at 

different frequencies to make a correction for a third. An even better 

situation is to have three surveys to fit the two parameters in Equation 

4.10. This would give an extra degree of freedom with which to test 

independently for a good fit. There are three surveys in the literature 

at 408, 6l0 and 1407 MHz with which to attempt to determine the rotation 

49 54 55 

measure in my region of interest. ? 3 

Although the resolution of all three surveys is approximately the 
same, 2° of arc, they don’t cover the same points. Because of this 
and because I have to wind up with a result which is averaged over my 
beamwidth, I decided to convolve the surveys with my antenna pattern 
before attempting to fit Equation 4.10. By doing so I felt that any 
systematics in individual readings or sparse coverage in a region for 
a particular survey would be averaged out in the large. Table 4.4 
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TABLE 4.4 

BINNED SURVEYS AT LOWER FREQUENCY 
BINNED 408 MHz SURVEY 



Bin 


Q 


u 


T 


°T 


Angle 


AAngle 


1 


0.2551 


-0.8107 


0.8499 


0.0206 


143.7336 


O.6947 


2 


1.2417 


-0.0071 


1.2417 


0.0215 


179.8353 


0.4954 


3 


0.4868 


0.2121 


0.5310 


0.0210 


11.7738 


1.1343 


4 


0.3362 


-0.l6l7 


0.3731 


0.0170' 


167.1560 


1.3049 


5 


-0.4497 


-0.2569 


0.5179 


0.0145 


104.8694 


0.8029 


6 


-0.7188 


-0.4034 


0.8242 


0.0148 


104.6497 


0.5137 


7 


-0.7658 


-0.2550 


0 . 8072 


0.0173 


99.2100 


0.6157 


8 


-0.0240 


-0.5514 


0.5520 


0.0223 


133.7520 


1.1585 


9 


0.1274 


-0.6139 


0.6270 


0.04'4i 


140.8636 


2.0131 


10 


0.3056 


-0.5857 


0 . 6606 


0.0330 


148.7786 


1.4292 


11 


0.5183 


-0.4465 


0.6841 


0.0221 


159.6279 


0.9245 


12 


O.1812 


-0.2787 


0.3325 


O.OI96 


151.5155 


1.6904 


13 


-0.4i66 


-0.0449 


0 . 4190 


0.0312 


93.0771 


2.1354 


14 


0.0855 


0.6023 


0.6083 


0.0159 


40.9592 


0.7475 


15 


-0.1455 


0.5694 


0.5877 


0.0115 


52.1664 


0.5588 


l6 


-0.1987 


0.5723 


0.6058 


0.0121 


54.5732 


0.5702 


17 


0.2000 


0.6266 


0.6577 


0.0131 


36.1509 


0.5706 


18 


0.3639 


0.5715 


0.6775 


0.0148 


28.7552 


O.6268 


19 


0.2257 


0.5482 


0.5928 


0.0160 


33.8108 


0.7746 


20 


-O.1929 


0.3432 


0.3937 


0.0151 


59.6691 


1.0996 


21 


-0.2014 


0.1726 


0.2652 


0.0176 


69.6953 


1.8996 


22 


0.1154 


0.0236 


0.1178 


0.0230 


5.7685 


5.5845 


23 


0.0574 


-0.0345 


0.0670 


0.0167 


164.4988 


7.1199 


24 


-0.3654 


0.1052 


0.3802 


0.0186 


81.9654 


1.4010 


BINNED 610 MHz 


SURVEY 










Bin 


Q 


U 


T 


ct t 


Angle 


AAngle 


1 


0.1924 


-0.0906 


0 . 2127 


0.006l 


167.3881 


0.8248 


2 


0.2718 


-0.0236 


0 . 2728 


O.OO67 


177.5216 


0.7058 


3 


0.2162 


0.0383 


0 . 2196 


0.0057 


5.0284 


0.7452 


4 


0.1281 


0.1774 


0.2189 


0.0047 


27.0834 


0.6l44 


5 


-0.0933 


0.1882 


0.2101 


0.0039 


58.1868 


0.5368 


6 


-0.2577 


-0.0665 


0.2661 


0.0039 


97.2307 


0.4185 


7 


-0.2955 


-0.0299 


0.2971 


0.0043 


92 . 8929 


o.4i4i 


8 


-0.1610 


-0.1064 


0.1930 


0.0050 


IO6.7299 


0.7448 


9 


-0.0889 


-0.1318 


0.1590 


0.0057 


117. 99^9 


1.0280 


10 


-0.7657 


-0.1950 


0.7901 


0.0057 


97.1432 


0.2085 


n 


-0.3553 


-0.2296 


0.4231 


0.0043 


106 . 4366 


0.2926 


12 


-O.1789 


-0.3978 


0.4362 


0.0034 


122.8918 


0.2224 


13 


-0.044l 


-0.1889 


0.1940 


0.0032 


128 . 4281 


0.4664 


l4 


-0.0413 


0.1135 


0.1208 


0.0031 


54.9980 


0.7461 


15 


-0.0885 


0.1691 


0.1908 


0.0032 


58.8083 


0.4837 


l6 


-0.1267 


0.1612 


0.2051 


0.0034 


64.0847 


0.4763 


17 


-0.5166 


0.1374 


0.5345 


0.0037 


82.5539 


0.1973 


18 


-0.6184 


0.3182 


0.6954 


o.oo4i 


76.3853 


0.1702 
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TABLE 4.4 (continued) 



BINNED 6lO MHz SURVEY (continued) 



Bin 


Q, 


u 


T 


°T 


Angle 


AAngle 


19 


-0.7566 


0 . 6826 


1.0190 


o.oo46 


68.9721 


0.1293 


20 


-0.4072 


0.4526 


0.6088 


0.0043 


65.9895 


0.2016 


21 


-0.1798 


0.2445 


0.3035 


o.oo4i 


63.1695 


0.3915 


22 


0.0294 


0.5464 


0.5472 


0.0042 


43.4574 


0.2201 


23 


0.1441 


0.3553 


0.3835 


0.0046 


33.9612 


0.3403 


24 


0.059^ 


0.0193 


0.0625 


0.0054 


8.9797 


2.4979 


BINNED 1407 MHz 


SURVEY 










Bin 


Q 


u 


T 


°T 


Angle 


AAngle 


1 


0.0524 


-0.0280 


0.0594 


0.0058 


165.9682 


2.8037 


2 


0.1347 


0.0405 


0.1407 


0.0058 


8.3572 


1.1818 


3 


O.ll4o 


0.0902 


o.i454 


0.0058 


19.1729 


1.1340 


4 


0.0185 


0.0619 


0.0646 


0.0055 


36.6985 


2.4473 


5 


-0.0260 


0.1522 


0.1544 


0.0056 


49.8548 


1.0467 


6 


-0.0926 


0.1388 


0.1668 


0.0055 


6l . 8627 


0.9396 


7 


-0.1452 


0.0527 


0.1544 


0.0053 


80.0338 


0.9861 


8 


-0.1256 


0.0104 


0.1260 


0.0055 


87.6273 


1.2446 


9 


-0.0540 


-0.0008 


0.0540 


0.0054 


90.4255 


2.8782 


10 


-0.0090 


0.0070 


0.0114 


0.0053 


70. 99^3 


13.3585 


11 


0.0195 


-0.0063 


0.0205 


0.0053 


170.9830 


7.3773 


12 


0.0054 


-0.0365 


0.0368 


0.0052 


139.2055 


4.0506 


13 


0.0113 


-0.0537 


0.0549 


0.0049 


140.9456 


2.5499 


14 


0.0843 


-0.0 4io 


0.0937 


0.0047 


167.0171 


1.4336 


15 


0.1101 


-0.0382 


O.1166 


o.oo44 


170 . 4431 


1.0888 


16 


0.0875 


0.0158 


0.0889 


0.0044 


■ 5.1105 


1.4255 


17 


-0.0199 


0.0815 


0.0839 


0.0044 


51.8526 


1.5014 


18 


-0.0685 


0.0232 


0.0724 


o.oo44 


80.6350 


1.7262 


19 


-0.0701 


0.0202 


0.0729 


o.oo46 


81.9476 


1.7981 


20 


-0.0522 


-0.0112 


0.0534 


o.oo48 


96.0603 


2.5859 


21 


-0.0465 


-0.0694 


0.0835 


0.0051 


118.0791 


1.7480 


22 


-0.0303 


-0.0769 


0.0826 


0.0053 


124.2484 


1.8538 


23 


0.0190 


-0.0859 


0.0880 


0.0050 


141.2228 


I.6261 


24 


0.0353 


-0.0920 


0.0986 


0.0051 


145.4948 


1.4855 
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gives the results of the convolution, where the convolved data have been 
binned in exactly the same way as my actual data. 

In attempting to fit the angles of these data to Equation 4.10, I 
fixed one of the three angles and tried the fit with various multiples 
of 180 degrees added to the other two angles . In this way I sought to 
account for possible rotations of greater than l 80 degrees between 

frequencies. I then selected the best fit f^or each bin based on the 

2 2 
X and have displayed the results in Table 4.5. In most cases the x 

is huge, though it should be less than one for a reasonable fit. An 

explanation is that the coverage in each individual survey is not 

good enough to allow accurate averaging. Coverage is an important 

factor in that the electron concentrations and fields change significantly 

on a small scale, at times even less than that of the resolution of the 

individual surveys. If not enough of the contributions are added, the 

average becomes inaccurate. Another possibility is that localized 

Faraday depolarization causes different areas at different frequencies 

to be important in the average of a particular bin. 

Because of the poor quality of the fits obtained for the three 
lower frequencies , I decided that the angles obtained from them would 
be useless for the purpose of correcting my survey and that I should 
forget them entirely and concentrate on working with the amplitudes. 

I felt that if there were a reasonable correlation between the amplitudes 
from my survey and the survey at the next lowest frequency, I could 
fit the amplitudes, using the spectral index as a parameter. 

Amplitudes obey Rayleigh rather than Gaussian statistics. The 
correlation coefficient and its associated expectation values for 
Rayleigh statistics are given by 
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TABLE 4.5 

RESULTS OF FITTING A ROTATION MEASURE TO LOWER FREQUENCY SURVEYS 



Bin "m ae m ’ x2 



1 


-32.82 


0.003 


14.16 


2 


63.18 


0.001 


1.92 


3 


63.22 


0.003 


21.44 


4 


-33.73 


0.006 


17.39 


5 


65.51 


0.002 


69.07 


6 


52.62 


0.001 
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<C ( 1 , 2 ) > = 0 

2 i 

<C(l,2) = W 



(4.17) 



For the correlation between the 3.2 and 21 cm data, I find that 
C(9370,l407) = .1668. Comparing this with tjie expected error for 24 
points shows that this is only about a .8 sigma effect and not a very 
convincing correlation. 

In spite of the low correlation, I still carried out the fit, 
feeling that though I could not use it for purposes of subtraction, I 
could at least ask whether or not the scale difference between the two 
sets of measurements was reasonable. The value of the spectral index 
which minimizes the Chi-square is -2.49, which doesn't seem too un- 
reasonable considering the possibility of Faraday depolarization. This 
confirms my back of the envelope calculation and indicates that on the 
basis of the 21 cm data I should be seeing an effect at 3.2 cm. To 
find the source of the lack of correlation, I plotted the scaled 
amplitudes of the two sets of data in Figure 4.5. The biggest 
discrepancy is the lack of a. significant contribution at 21 cm between 
19 and 23 hours of right ascension. In this area is the Cygnus arm of 
the galaxy, the hottest region that falls in the beam. 

Going back to Table 4.4, there doesn't appear to be a contribution 
in this region at 4o8 MHz either, yet looking at a survey at the same 
frequency and an 8° bearawidth, there is an indication that an effect in 
fact does exist. In this region especially, where the galactic 
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magnetic fields are complex and are directed along the line of sight, 
one would expect Faraday depolarization and problems of coverage to be 
important . 

I feel that there is strong evidence for a contribution to the 
3.2 cm polarization from the galaxy. Unfortunately, the data available 
at other frequencies does not seem capable of supplying corrections. 

In all fairness I can understand why they wouldn't, given the complex 

i 

nature of the emission regions and magnetic fields in the galaxy, the 
large gap in wavelength between them and 3.2 cm , and the purpose for 
which these surveys was made, which was to provide high resolution 
coverage of the most important regions of polarized radiation at their 
frequencies. In their case the most important regions were close to 
the galactic center and therefore their coverage was concentrated in 



this area. 
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4.5 FINAL ANALYSIS 

Being unable to subtract the galaxy, I have to accept the harmonic 
components as they are given in Table 4.2 and investigate the upper 
limit they infer. From here the problem lies in the complicated 
predictions of the model given by equation 2.l6. Rewriting this 
equation twice, once with (3 = 0 and once with (3=45°, I produce the 
predictions for the Q and U parameters, a*t a declination 6 

Q = (T - T ) [sin2§sin6cos6sin(t-0' -rr ) + sin §(l-icos 6)sin2(t-Q' +tt )], 
w a max o o -jj- 

2 

U = (T - T ) [sin2§cos6sin(t-a +tt) + sin §sin6sin2(t-0' )] , 
v w a'max v o ' v o' ’ 

(4.18) 

where everything is expressed in terms of sine functions. These may be 
compared directly with my data. The problem in doing this lies in the 
fact that the phases of the experimentally measured components don’t 
agree well with the model. This isn’t surprising, considering that most 
of the contribution is due to galaxy and not an asymmetrically expanding 
universe. For example, the measured phases of the 24 hour components 
for the Q and U parameters, when compared with Equation 4.13* predict 
values for ot which differ by 150.24 degrees, well outside the exper- 
imental errors . 

This suggests that by forcing a fit to the model, the non-agreement 
of phases would help toward a lower limit. I object to this course for 
two reasons. First the error introduced by the phases may be large 
enough that any lowering of the limit would probably be offset by an 
increase in the uncertainty. More importantly, I feel that because I 
do have statistically significant amplitudes, one can always argue that 
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even were they to be predominately from the galaxy, they might be 

masking a significant cosmological contribution. A fairer estimate is 

to compute a 9 0 % confidence upper limit for each 12 and 24 hour harmonic 

of Q, and U and forgetting phases, compute a result in each case for 

the absolute value of (T - T ) . This procedure assumes that the 

v w a max 

phases in the data fit the model perfectly and uses the amplitudes 
to produce an upper limit. Another advantage of this method is that 
the statistics of the amplitudes are well known and setting confidence 
levels is relatively easy. First we form the quantity, 



where T is assumed to consist of a random component plus a real signal, 
and the expectation value is based solely on random errors. This is 
just the signal power divided by the noise power, a parameter which 
has a non-central Chi-square distribution for two degrees of freedom. 

64 

Consulting the appropriate probability graph for this distribution, 

I find the following upper limits 




(4.19) 



24 hr: T < .75 mdegK 



For Q, 



For U 



12 hr: T < .38 mdegK, 

24 hr: T < .8l mdegK 

12 hr: T < .67 mdegK, with 90 io confidence. 



(4.20) 



Inserting these into the expressions for Q and U along with the 



declination of the observation, 40.35 degrees, one obtains 
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For Q: I (T - T ) I < 1.6l mdegK, 

1 v w a ; ma x 1 5 

U: | (T - T ) | < 1.31 mdegK, with 90 % confidence, 

w a max 

We must also consider the possibility of $ being very small. In 
this case the cosmological contribution would be buried in the DC 
contribution. Assuming ® = 0 and using the upper limit on the DC 
level, I find from Equation 4.1 




|(T - T ) | < 1.47 mdegK, with $ 0 % confidence. (4.22) 

w a max 

Taking the largest of the three, I find for an upper limit on the mag- 
nitude of the polarization of the microwave background: 

< 1.6l mdegK, with $ 0 % confidence, (4.23) 

In terms of a 2.7 degree microwave background, this amounts to . 06 % 
polarization. Inserting this result into Equations I.19 we obtain for 
e the asymmetry parameter and A h o /h Q : 

|e | < 8.15 x 10 _1+ 



A h 



- < 6.4i X io“ 5 . 



(4.24) 



h ~ 
o 



_ rr o 

where I have assumed the nominal values n Q = 10 /cnr and h^ = 50km/see 
Mpc \ It should he remembered that Equations 1.19 were derived assuming 
that l(t), the ionized fraction, was equal to one for all times. In 
general this is an unreasonable assumption, but Rees has shown that in 
the event of a reheating phase at a redshift of z > 7 (due perhaps to 



galaxy formation), that the solution applies and defines the upper limit 

42 

on the expansion asymmetry. 

In the event that there was no reheating between recombination at 
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z 1000 and the present, the power anisotropy measurements imply a 
polarization smaller than that given in 4.l8, and in this case 
provide the upper limit. 

4.6 FUTURE DIRECTIONS 

As in any enterprise, there is always room for improvement and 
this one is no exception. The two most glaring limits on this exper- 
iment are the galactic contribution and the effects of ground radiation. 
In principle, the latter could be solved by improving the ground shield 
and by using an antenna with somewhat better sidelobe characteristics. 
In the case of the galactic contribution, there doesn't seem to be any 
reasonably clean way of correcting for it at 3.2 cm. For this reason, 

I feel the best solution is to try a frequency where the expected 
contribution from the galaxy is less than the expected limit from 
receiver noise. For example, at 1.6 cm and assuming a spectral index 
of -2.5, the galaxy emission should be down by a factor of 5*7. If 
one were to use a receiver of comparable noise figure to the one used 
in this experiment, the galaxy would be well down in the noise, while 
the power from the microwave background, having a blackbody spectrum, 
would be roughly constant. If there were still no evidence of a 
cosmological polarization, the improvement in the limit would be 
roughly a factor of two. 

Another major improvement oould be made by using the second part 
of the dual mode transducer on the switch in a two receiver system. 

This would be a form of Graham's receiver which would decrease the 
effective system noise by the square root of two. Because the signal 
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in the two arms of the transducer are out of phase at the switching 
frequency, while the switch offset is in phase, it also offers the 
opportunity of balancing out the switch offset before entering the 
lock-in. In this way the immunity to gain drifts is enhanced. 

Taken together, these changes would give ah improvement of 2.5 
in overall sensitivity to the background polarization. 
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APPENDIX A 



Al. TIME DEPENDENCE OF THE SCALE FACTORS IN AN AXIALLY SYMMETRIC 
EUCLIDEAN UNIVERSE WITH NO MAGNETIC FIELD 
In order to obtain the time dependence of the scale factors, one 
has to obtain the solutions to Einstein's equations for the metric in 
question, given an appropriate stress tensor^. In this work I deal 
with an axially symmetric Euclidean metric given by 
ds 2 = dt 2 - A(t) 2 (dx 2 + dy 2 ) - W(t) 2 dz 2 . Metrics of this type 

40 

have been considered by many authors, including Thorne (1967), D 
who has given the following results and whose notation I follow. 

In the absence of magnetic fields, I assume a stress energy 
tensor of the form, 



/ =rp 0 0 O-i 

V 

0 -p/c 2 0 0 

0 0 -p/c 2 0 

u 0 0 0 

for use in the equations. The equations for a quite 

57 

have been written down in detail by Dingle (1933). 
for the above metric we obtain 



-p/c 2j 



(Al.l) 



general metric 
In particular , 



8ttGP _ A W AW 

4 _ A + W + AW 5 
c 




8rrG 



•• • p 

2A A _ 

A + 2 ’ 

h k 



2m 

A 2 + AW 



(2 equations) 



(A1.2) 



89 



By subtracting the first equation from the second we find the following 
expression. 



• • •• • p • • 

AW A_ AW _ 

A " W + A 2 " AW 



(A1.3) 



for which W = A is clearly a solution. In order to find others, we will 

make the following two assumptions. Since in the limit of A = W, we have 

an Einstein-de Sitter uni verse for which A - A^t/t^) 2 ^, and since the 

case we are considering is a small simple departure from the symmetric 

case , I will assume that A is a power law in t, i.e. A = 

In addition, I will assume that W is of the form W = A + qA^, where 

q is some arbitrary constant. In effect this says that for small 

asymmetry, we expect a solution. close to that of the Einstein-de Sitter 

universe and for that reason I choose to expand A in terms of the 

Einstein-cfe Sitter solution. The choice of q for a constant is 

unfortunate and shouldn't be confused with q^, the deceleration 

parameter. I use it to provide continuity with Thorne's results. 

, N 

Inserting A + qA for W into Equation A1.3 we find 

q[(l-N) A W_1 A + (1-N 2 ) A N-2 A 2 ] = 0. > (A1.4) 

Inserting the explicit form for A gives 




r / _ \ , / .. TT \ 2, aN-2-i _ 

La(a-l)t + (1+N)cv *c J = 0, 



(A1.5) 



which reduces to 



(N+2)a = 1. 



(ai.6) 
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Next, taking the third of Equations A1.2, we note that the left 
side depends on the density, which is inversely proportional to the 
volume, meaning that aSj is a constant. Multiplying this equation 
through by A 2 ^ thus gives 



A Sf + 2AAW = constant, 



(A1.7) 



which becomes, after substituting for W: 



3AA^ + qA^A^ (l+2N) = constant. 



(A1.8) 



Again substituting for A and using the result A1.6, we obtain the 
expression for a, 



^3a-2 /3a-2 N / o \ - -3 a -2 

3t - (— — ) q ( — ) a t = constant, 



(A1.9) 



which is constant if 3<*-2 = 0 or a = 2/3. This, with equation A1.6, 
implies that N = l/2, giving the following forms for A and W: 



A = A (t/t ) 
o w o' 



2/3 



W = A + qA 



- 1/2 



(A1.10) 



From the above we can now produce an explicit form for the 

differential Hubble r s constant: 

W A 

to = W - A • 

By plugging in the appropriate expressions for A and W, we find 

3/2 



<iK 



Ah = 



3t 



A Sj 



(Al.ll) 



(A1.12) 



which gives us the result that A WAh = Ah^ is a constant. 



A2. THE ANGULAR DISTRIBUTION OF THE TEMPERATURE OF RADIATION IN AN 
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ANISOTROPIC EUCLIDEAN COSMOLOGICAL MODEL 



This appendix will demonstrate 
anisotrcpically expanding Euclidean 
in the temperature of the microwave 
we mean models with a metric of the 
W(t) 2 dz 2 . 



how, in the absence of interactions, 
models give rise to an anisotropy 
background'. By Euclidean models 
form ds 2 = dt 2 - A(t) 2 dx 2 - B(t) 2 dy 2 



In order to effect this, it is convenient to first show that 

p, = U, = dx,/dT is conserved over the trajectory of each individual 

blackbody photon. In the standard way, I define the metric tensor 
i v 

as g , in terms of which the affine connection is given by 



^ = l/2q°^ [g , + g , - g , ] , 

Py Pp y py P yP p 



(A2.l) 



in terms of which the geodesic equations of motion of a photon may be 
expressed as 



du j. O’ B 

-r— = - 1 . U U r . 

dT ap 



(A2.2) 



From these definitions 3 the evolution of along the trajectory 
of a photon may be obtained. Consider first 



dT 



du. 

1 

dT 

d 



dt <%Y* V > 



Yd, s du 

“ U dT ®iy + e i Y dT 



Y 



(A2.3) 
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Substituting from Equations A2.2 we get 



dP. 






(A2.4) 



i Y B / _Y ( 

dT ^.Y ’ p u u + ®Ly r ap U 

which by the definition of the affine connection, Equation A2.1 becomes 



dP. 



a 6 

n r 



6> p + %p> Y ' V " 



dT = > p uV ‘^Y (gY6) 

- hf • p aY " P ‘ i %i > p u “ u? - i %p > « u “ u 4 * > i u<,ulS - < A2 - 5 > 

In the above, y is only a dummy index and may be replaced with a , giving 



dP. 

dT 



i / i yOBi dBi a B 

r = ( %c > p - 2 . p ) " " - s 

y* or B 
= r . u u . 



aiB 

Because dp^/dT has to be symmetric i-n a and B > which are just 
dummy indices, it must be equal to the same expression as follows: 



(A2.6) 



dP. 



dT 



- = Hr . Q + r ] u 



a B 



Q'iB Bi a 



u' 



i ri 



Cig,^ o + s^, , - gift*) + i (g. Q , + g . . - g. , )] 



a B 



iOf B B a ' i 
i , x a b 

■ “ ( V i + V i> u “ 
i of B 

= * V i uu • 



iB Of ap i ia B 



(A2.7) 



lu u 



The last line is the needed result. It states that if the metric 
tensor is not a function of x\ then = I is a conserved quantity 
along the trajectory of a photon 5 in the absence of scattering. In this 
work I deal with metric tensors which are functions of coordinate time 
alone , which implies that the photon momenta p^> i = 1, 2, 3? are 
constants of the motion. This key result allows us to compute the 
redshift of a photon emitted at an epoch t and observed at the present 
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time t from the direction 'n\ Starting with the original energy of the 
photon, as measured by an observer at rest ’.rith respect to the x, y, and 
z coordinates, we can write 
(hv) 2 = E 2 



= P 

= P g X ^P + P g y ^ P + P g 2 ^ P 
x s X y 6 A z s V A 



= - hr"- (p x ) 2 + ~r~ (rJ 2 + -j— 0? z ) 2 ] , 

X B 2 (t) y C 2 (t) Z 



A 2 (t) 



(A2.8) 



which by conservation of the momenta p^ is equal to 



(hv) 2 = - [~ 2 ±~ (P 0 ) 2 + —^~ (P 0 ) 2 + -| (P 0 ) 2 ]. 

A 2 (t) x B 2 (t) y C 2 (t) 2 



(A2.9) 



For each component of the momentum we can show the following : 

? A 



(V 2 - P .o 



X 



= P o P 
x u 

- p 2 (fr • / xj 2 

o v ' 

- (hv,) (-S' • S') 2 , 

which taken together with Equation A2.8 gives 

,2 



(hv) = - 



(A2.10) 



B 2 C 2 

('n' • / x) 2 + -2 fn' • ^) 2 + -§ — ('fr • 'x') 2 

A (t) B 2 (t) 



c 2 (t) 



(hv ) 2 . 
' o 



By rearranging terms and taking square roots we arrive at the following 

58 

equation for the redshift z: 



9 ^ 



1+z v 

2 2 - — 

. 2 V 2 2 2 

( / n v • / x) + ("fr • 'y') + — ('n' • / z') 

B (t) C 2 (t) 

(A2.12) 

With this result, the temperature of the microwave background as a 

function of the direction of observation can finally be computed. The 

argument used is similar to that given for tjie isotropic case by 

Tolman and rests on the proof of Liouville's theorem for general 
59 

relativity. ' Liouville's theorem states that in the absence of 
interactions, the density of states, N, along the geodesic of a particle, 
is a constant. Assuming that we are considering a universe containing 
radiation with a blackbody distribution at recombination, or at least 
at the epoch of the last scattering, implies 




N = 



a hv/kT 

2 ~ J- 



e 

= N 



hVo/kTo.! 



(A2.13) 



using the Planck distribution characteristic of a blackbody. 
Equality here implies equality of the arguments of the exponentials 
giving 



v 

V o 
T ~ T 



(A2.14) 



or 



v T 
o o 



1+z 



V 



(A2.15) 
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Combining this with Equation A2.ll, one obtains the final 
expression for the angular distribution of the temperature of the 
microwave background in an anisotropically expanding Euclidean universe, 



B 



T = 7^— = T [-2- sin 2 0 cos 2 0 + sin 2 0 sin 2 0 + — 2- 






1+2 



COS 0] 2 . 



B 



(A2.l6) 

The result for an axially symmetric universetis a specialization with 

A = B and C = W 2 

A ? W i 

T = T [-2- sin 9 + -2- cos 0] 2 , (A2.17) 

° A W 



manipulation of which gives 



W r 

T = T £ [1 + 

o W . 

o A 



A 2 !-, 2 



2 H 2 



- 1 sin 2 0] 2 . 



(A2.18) 



In the event that the anisotropy in the expansion is small or, in other 
words , for A 2 W 2 /A 2 t'7 2 = 1, the binomial theorem allows us to write 



p 

T = T (l + e sin 9) , 
o av v ’ 5 



(A2.19) 

where T = TV/ /w and e = h (l - A 2 W 2 /A 2 V-/ 2 ), a small number . 
av 1 o * s o ' o 

From Appendix A1 we can write the explicit expressions for A and 



W 



A = A (t/t ) 2 / 3 
o ' ' o 

l 



(A2.20) 



W = A + q A" 2 . 

Insertion of these into Equation A2.17 allows us to compare the extent 
of the expansion asymmetry with the asymmetry in the temperature 
distribution, 
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e = 



2A~ 



q (1 - ( t o/t) 3 ) + 2 A q 3 / 2 (1 - ( t o/t) 3 / 2 ) 

W 2 R 

1 + 2q/A + q /A ° 



. (A2.21) 



o - 1 ' o 

From this it is clear that the requirement that e be small, so that 
the binomial expansion in Equation A2.l6 is valid, is equivalent to 
the parameter q being a small number and thus, the asymmetry in the 



expansion being small, which is exactly our above assumption. 



A3. EFFECT OF A SINGLE THOMSON SCATTERING ON THE POLARIZATION OF THE 
MICROWAVE BACKGROUND IN AN AXIALLY SYMMETRIC UNIVERSE 
When the energy of a photon is much less than the rest energy of 
an electron, its primary interaction with ionized matter is through 
the mechanism of Thomson scattering. This elastic form of scattering 
can be looked on as radiation by an electron which has been accelerated 
by an incident electromagnetic wave. Since this is electric dipole 
type radiation, the polarization of the outgoing wave can be derived 
from the solution for a simple radiating dipole. If the acceleration 
of the electron were considered to be in the a direction, the polar- 
ization of the outgoing wave would be perpendicular to the direction 
of propagation and would lie in the plane containing the direction of 
propagation and the direction of the acceleration of the charge (see 
Figure A3.l). To obtain Ie', the unit vector in this direction, we 
make the following definitions: 

"r s = (sin9^, 0, cos9^) , the line of sight to the observer 
'a = (sin9 cos sin9sin$, cos9), the direction of the 
acceleration of the charge. 

t = the angle between / r v and / a\ 



m> 
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DIRECTION OF SCATTERED RADIATION FOR 
THOI4SON SCATTERING 
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From the description given above of the constraints imposed on 
the polarisation of the scattered radiation, we can write down three 
equations for the three components of "S'. Perpendicularity to the direction 
of propagation / r > requires that ^ • 'r' = 0, while lying in the plane con- 
taining / r' and 'a' means that Is' • ( / r v x ‘si) = 0. Finally, being of unit 
2 2 2 

length implies E^ + E^. + E z =1. 

Writing these out in detail, we obtain \ 



£ • * 



r = E sin9. + E cos9. = 0, 
x 1 z 1 



E * (r x a) = -cos9, sin9 sin§ E + (sinC, cos9-cos9 sin9 cos?) E 
' 1 x ' 1 J 



+ sin9, sin9 sin§ E , 
1 z 

2 2 2 
E+E + E = 1. 
x y z 



(A3.1) 



Another expression which will prove useful is 

sin^ty = |r x a |^ = sin^9 sin^$ + (sin9^ cos9 - cos9^ sin9 cos§)^. 



Solution of the above equations yields the following for the 

components of E 
2 

_ cos 9 (sin9 cos9 - cos9 sin9 cos§) 

E ^ = , 

X .2 

sin 



(A3. 2) 



2 2 2 2 
Ey = sin 8 sin §/sin \jf, 

2 

sin 9 (sin9 cos9 - cos9 sin9 cos§) 

E = — 

z . 2 . 

sin ijf 



(A3. 3) 



which can now be used to describe the polarization of the scattered w ave 
In order to put this information in a more useful form, I shall follow 
Rees' notation and define the polarization of the outgoing w ave in terms 
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of the plane containing the direction of propagation, r and the symmetry 

A 

axis of the universe, W. The direction perpendicular to this plane 
will be called the 'a' direction, while vectors lying in this plane 
and perpendicular to r will be said to be in the 'w' or cqplanar 
direction. Applying these definitions to Equations A3. 3, we find 



E = E 



y 

2 2 2 
sin 9sin §/sin i|r , 



E 



w 



2 2 
E + E 
x z 



2 2 

= (sin0^ cos0 - cos0^ sin0 cos§) /sin ijr 



(A3. 4) 



If the incident radiation is polarized in the T a* direction and 

has direction of incidence 0 , the acceleration of the charge will 

lie in the x, y plane and be perpendicular to the direction of 

incidence. Thus 0 = tt/2 and $ = 0^ ± tt/2. If it is in the I w’ 

direction 0 = 0^ ± tt/2 and $ = + tt. Taking these cases separately , 

the two components of scattered radiation (as seen by an observer 

whose line of sight r, makes an angle 0, with respect to the axis of 
A n 

the universe W) will be computed. 

In Appendix A2 it was shown that in an anisotropically expanding 

axial universe with a small anisotropy, the temperature of the 

radiation can be characterized as having a distribution of the form 

2 

T = T^ v (l + e sin 0), where e is an asymmetry parameter and 0 is 
the angle between the direction of observation and the axis W. In 
this case, in order to deal with the polarization, define two 
asymmetry parameters and and specify the temperature of the 
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radiation in a given polarization in terms of them by writing 
2 2 

T = T(l + e sin 0 1 ) and T = T (l + e sin 9.. ) . Any change in the 

cl cl X W W _L 

polarization due to scattering or other effects can now be characterized 

by its effect on these parameters. 

The amount of radiation scattered in any given direction is given 

cLcr 3 2 

by the product of the incident flux and — = ^ o^sin K the 

differential Thomson scattering cross section, where o^ is the total 

scattering cross section and \|f is the angle between the direction of 

acceleration of the charge and the direction of scattering. A fraction, 
2 

E , of this is scattered with polarization in the ! a r direction and a 
a 

2 

fraction, E^ , with polarization in the 'w' direction. Taking the 'a' 
direction first, the amount scattered into this polarization from 
incident radiation, initially polarized in the same direction is 
given by 



T ' = P T ~ E 2 dO 
a_^ *) a dfi a 2 

■?T o 2 2 

D T p /, . 2 n s . 2, sin 9 sin 5 _ n 

J (1 + e sin 9 2 ) sin iji d cos0^ d<J>^. 



8n 

3T ct_ 



- 2 i 

sin v 



2 2 



= — 7 s — — P (l + e sin 2 0 o ) cos^$ 0 2d cos0 o d$_ 

OTTtJ N a 2 2 22 



(A3. 5) 
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In a similar way, radiation scattered into the 'w' polarization is 



given by 
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T ' = f T ^ S 2 dO, 
w J a dQ w 2 



T 

8n 


J ( 1 + 


e sin' 
a 


3T a T 
8 tt 


J ( 1 + 


. 2 
e sin 1 
a 


Tct t 


2„ 


, 3 
( 2 + 


2 


cos 0^ 



2 (sin0^cos0 - cos0 sin0 cos§Y 

1 * 

sin ^ 



d cos 0 2 d $ 2 



a' 



(A3. 6) 



Starting with incident radiation polarized in the 'w' direction, 
the corresponding results are 



T = f T E 2 dQ 
a„ J w dQ a ; 



w dQ a 2 



3T a T 
8rr 


J a- 


3T a T 
8rr 


I a - 


OJ 

EH 

11 


(1 - 


for radiation scattered 


polarization, 




T ' = f T ^E 2 , 

J v dQ *w 


3T a T 
8 tt 


Ja 


3T a T 
8n 


ja 


11 

i -3 

ro 


HlCVJ 



■ 2 , 



2 2 

2, sin 0 sin $ 
i ♦ 

• 2 , 

sin \|r 



d cos0 2 d$ 2 



2 2 2 
s^sin © 2 ) cos ©2 sin 

f e ) , 

5 ir ’ 



(A3. 7) 
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0 _ (sin0 cos0 - cos0 sin0cos$)‘ 

• /"V \ • U I J. 

;in ©2) sin y 



■ 2 . 

sin y 



dcos02d$2 



f (l + e sin 2 0 o ) (sin0. sin0^- cos0. cosO^osqO^d cos0 o d§ o 
J v w 2 /v 1 2 12 2' 22 



I %) + 1^1 si n 2 0 l( | ♦ 1 e„). 
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Adding these results gives the total contribution to each 
polarization 

T a' - \ + \ = To T f 1 + i V io £ w ) 

T w’ - T v’ + V, = To T C(1 OV 15 c „> + < 30 V i c a )sil, \ ] ' 

(A3. 9) 



If we now re-express equations A3. 9 in terms^ of new asymmetry parameters 

e 1 and e ' we obtain to first order in e and e ,the following sol\itions 
aw a w’ e 

for the new parameters : 
e ' = 0 



(A3. 10) 



e ' = — e - i e , 
w 30 w 2 a ’ 

where the new average temperature is T' = Ta^,(l + \ e & + i_ e w ) • 
In matrix form this can be written 




(A3. 11) 



This last expression is the final result for the effect of a single 



Thomson scattering. 
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a4. terrestrial observation of a cosmological linear polarization 

In previous sections I have shown how cosmological causes can give 

rise to a linear polarization of the microwave background of the form 

(T - T ) = (T - T ) sin 2 0, (A4.l) 

v w a.' v w a max ’ v ' 

where 0 is the angle between the line of sight of the observer and 

the direction defined by the symmetry axis of the universe. Now, the 

question is what form does this polarization^ take for a ground observer 

trying to measure its effect? In my case I will assume the following 

measurement scheme. A microwave polarimeter will be aimed at the 

zenith and, as the earth rotates, will be allowed to sweep out a circle 

of constant declination in the sky. Given this situation, we need to 

find the variations in the signal produced by the rotation of the 

earth. Vie can diagram the experimental setup as in Figure A4.1, 

where E^ and E^, the unit vectors in the perpendicular and coplanar 

directions as defined in Appendix A3, have been drawn in. 

By stating their relationships to the directions R and W and to 

the normal to the R, W plane R x ¥, we can develop equations for the 

determination of E and E . First, because both of .them are parallel 

aw’ * 

to the polarization of incident radiation, they must be perpendicular 
to the line of sight, R. Next, E^ lies in the R, W plane and must be 
perpendicular to its normal while E is perpendicular to this plane 
and thus must be perpendicular to W. Finally each of them must be a 



unit vector. In vector notation these conditions become 
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Z N 




N- NORTH CELESTIAL POLE 
W- AXIS OF THE UNIVERSE 
R- LINE OF SIGHT OF THE OBSERVER 
8 = 7T/2 -S 

ot= cl 0 -t (where r is local sidereal time 

in angular units) 



FIG. A 4.1 

RELATION BETWEEN THE AXIS OF THE 
UNIVERSE, THE ROTATION AXIS OF THE 
EARTH AND THE LINE OF SIGHT TO THE 
OBSERVER 
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E • fi = 0 

W 



£ *8 = 0 
a 



E • (W x fi) = 0 
w 

* 

2 2 2 
E + E + E 
wx wy wz 



= 1 



g • W = 0 

a 

2 2 2 
E + E + E =1, 
ax ay az 



(A4.2) 



while written out in detail, they are 



E sin 6 + E cos 6=0 
wx wz 



E sin6 + E cos6 = 0 
ax az 



E sin $ sina cos6 + E (cos9sin6 
wx wy' 

-sin$coso'cos6) + E sin§ sina' sin6 = 0 
wz 

_ 2 2 2 
E + E + E =1 
wx wy wz 



E sin§ coscv + E sin9 sina 
ax ay 

+ E cos§ = 0 
az 

2 2 2 
E + E + E =1. 
ax ay az 



Taking these with the expression 



(A4.3) 



2 

| W x R | = sin 0 = sino sin a + (cos$ sin6 - sin§ cosa' cos6) 



2 - 






(A4.4) 



leads to the following solutions for the unit vectors; 



E 



cos6 (cos§ s in 6 - sin$ cosa cos6) 



vx 



sin9 



ax 



c os6 sin$ sincy 
sin9 



E = + 
wy 



sin$ sincy 
sin9 



E 



ay 



s in 6 cos$ - sin$ cos# cos 6 
sin9 



„ sin6 (cos£ sin6 - sinf coso' cos6) 

E = + r- 7 ' 

wz sin0 



E 



az 



sin6 sin9 sine* 
sin0 



(A4.5) 



this allows us to express the resultant field seen by the observer as 



E = E E + E 
a a w 



a 



E 

w 



(A4.6) 



with the amplitudes E^ and E^ complex numbers and functions of t. In 



order to resolve these into components relative to directions with 
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respect to the north celestial pole, we notice from Figure A4.1 that 
our choice of coordinates implies that the East-West component of 

— i 

E is the E component while the Worth-South component is given by the 

y 

projections in the plane of the other two or 



E = E sin 6 - E cos6 
ns z x 



E = E . 
ew y 



(A4.7) 



\ 



By applying a rotation to these two components , the field can he expressed 
relative to any direction we choose. For an axis at an angle (3 with 
respect to North > we can give the components with respect to the new 
axis as 





ns 



E 



(A4.8) 



ewi 



For a polarimeter of the type I am using, the important quantity 

2 2 

is , which is equal to c(|\l - I Eg | ). The quantity in brackets 

can be written in terms of E and E as 

ns ew 

| El | 2 - |E 2 | 2 = (I E ns | 2 - |E ew | 2 ) cos2p ♦ 2 ^ | sin2 Pj 

(a4.9) 



or in terms of E and E as 
a _ w 
■2 i- |2' 



IeJ 2 - |e 2 | 2 = 



( pl - K n 



[(cos§ sin6 - sin§ coscv cos6)^- sin^§ sin^cv] 



sin 0 



* 



2 |E -E , 

4 (cos§ sin6 - sin§ coscv cos6)sin§ sina 

sin 0 



cos2p 



42 



( l E I'- |e I 2 ) 

' 1 w 1 1 a 1 

2 

sin 0 



(cos§ sin6 - sin$ cosu cos6) sin$ sina 



E . E 

1 w a 1 r p p , _ __.p„ 

+ n Lsin § sin a - (cos§ sin6-sj.n§ coso- cos6) ] 

sin 0 



sin?f3. 



(A4.10) 
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In a practical situation where the signal is integrated for times much 
greater than the reciprocal of the receiver bandwidth, we are taking a 
time average of the above expression. Because the signals in the two 
polarizations, E & and E , are independent, the time average of their 
product is zero and we obtain 

Gfcos$ sin6 - sin§ cosa' cos6) - sin § sin a] cos2B 

\ (A4.ll) 

+ 2(cos§ sin<5 - sin§ cosa cos6) sin§ sina sin2p| , 

where we have multiplied through by the constant C and made the sub- 

2 

stitution (T - T )/sin 0 = (T - T ) . This result is the signal 

v a w'' v a w'max 

desired. Replacing a and 6 with their equivalents in terms of the 
local sidereal time t (this is equal to the right ascension of the 
zenith), the right ascension of the axis of the universe, a 1 , and 
the declination, 6 , and doing some manipulati on we arrive at the 
following: which is the polarization signal received by a switched 

polarimeter pointed at the zenith which has its positive axis making 
an angle (3 with respect to North. Here, positive is defined as 
the direction: 

T - T_ = (T - T ) )[cos^6 (l - •=! sin^a>) + sin^$ (l - ^ cos^6) cos2 (a - t) 

1 2 v w a max I K 2 ' v 2 v o 

1 2 

- — sin2§ sin26 cos^^- t)] cos2(3 - [sin^$> sinS sin2(o? o - t) 

- sin2§ cos6 sin (or - t)] sin2p 



(A4.12) 
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APPENDIX B 



Bl. THE ANTENNA 

As I had nothing to do with design of the antenna, I will not go 
into detail on the techniques of design but will only present a 
description of its operating characteristics. Interested readers can 

6 ( 

refer to the references , in particular the pa^er by Roll and Wilkinson, 
for further information. 

The antenna used was a conical, optimum gain horn with an aperture 
choke for suppression of the backlobes. Provision was made as shown 
in Figure Bl.l for the insertion of a calibration probe in the side of 
the horn, and this feature was in place during about two months of 
observations. The full width between half power points is about 15 
degrees which gives the polarimeter a resolution of one hour in right 
ascension. Figure B1.2 is a plot of antenna gain as a function of angle 
from the horn axis for illumination by orthogonal polarizations. As 
the procedure varies slightly from that used for a power radiometer, 
a little explanation is in order. In making the measurement, the 
polarimeter was mounted horizontally on a rotating table with the 
apex of the horn at the center of rotation. One polarization sensitive 
axis of the polarimeter was aligned horizontally and the other 
vertically. Incident radiation was provied by a klystron radiating 
through rectangular waveguide and a standard gain horn, which gave a 
linearly polarized source. Horn to horn distance was a little over 
twenty meters, giving a good approximation to a plane wave at 3.2 
centimeters. Two scans were made, one with the incident polarization 
horizontal, which corresponded to the E-plane measurement of a single 
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ANGLE FROM HORN AXIS 

FIG. B 1.2(a) 
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RESPONSE OF 3.2cm POLARMETER TO 
UNPOLARIZED RADIATION ( E PLANE - H PLANE ) 
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polarization microwave horn, and the other with polarization vertical, 
corresponding to the H plane measurement. Because of these analogies I 
have used the terms E and H plane in labeling the resultant plots. The 
appearance of the E plane measurement is very similar to what one would 
find for a power radiometer, while the H plane measurement has a striking 
difference. For some angles the sign of the signal changes, indicating 
a coupling between the polarizations. It is ^iot known whether the 
observed effect is due to external reflections or to coupling within 
the polarimeter itself. Because the effect occurs at levels of -35db 
and below, it has little effect on the experiment and was thus not 
pursued further. 

The primary purpose of the measurements was to produce Figure B1.3, 
which is the difference between the E and H gain patterns. This 
difference gives directly the response to unpolarized radiation, which 
can then be used to compute the expected contribution from ground 
radiation. The antenna temperature due to any source falling within 
the beam pattern is°^ 

T a = J* TP(Q)dn/Jp(G)dfi, • (Bl.l) 

where T is the equivalent blackbody temperature of the source at 3.2 
centimeters and P(Q) is the gain of the antenna. The denominator can 
be evaluated by noting that calibration is performed by covering the 
entire antenna pattern with a source at temperature T and calling the 
resultant antenna temperature T. For a rough approximation of ground 



radiation we can use the form 



11 4 



V TP Av( n ) An Af P ( n ) da 



(B1.2) 



At the observation site on the roof of Jadwin Physics building, 
the most prominent terrain feature and the expected source of most of 
the ground radiation is a structure on the roof, adjacent to the polarimeter. 
This structure houses the elevator motors and is approximately two 
stories high or an angle of about a half radian above the horizontal. 



Its width is about a third radian , giving it a size of about one 
sixth steradian located between an angle of sixty and ninety degrees 
from the horn axis. In this position the antenna gain from Figure B1.3 



as an appropriate contribution to the signal from ground radiation for 
an unshielded antenna. 

During the experiment the antenna was shielded from ground radiation 
by a shield which extended from below a zenith angle of ninety degrees 
to a zenith angle of thirty degrees. Measurements of its effectiveness 
using the klystron source were done for zenith angles of between eighty 
and ninety degrees. The average reduction in sensitivity produced by 
the shield in this range was 12.3db. If we take this figure and apply 
it to the above calculation, we find, assuming that the shielding factor 
is valid over this range, that a reasonable figure for the expected 
contamination from ground radiation is 7-9 mdeg. K. In the first two 
months of the experiment I noticed a DC level of polarization of 6.66 
mdeg. K., which agrees with the rough estimate above. I felt that this 



is about -37db for unpolarized radiation. If we take 290°K as an 



average temperature for the structure we obtain 




TP poll (0)Afi 




(B1.3) 
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was probably not unreasonable. An order of magnitude measurement, using 
the shielding factor and a piece of eccosorb as a source inside the 
shield, showed rough agreement. In fact, this measurement was correct 
at that time. It wasn't until I later decided that I would not need 
the calibration probe and therefore removed it from the horn that I 
discovered that the high degree of sensitivity to ground radiation was 
due to the perturbing influence of the probe . s In the months that 
followed the level of DC polarization dropped to less than 1.5 mdeg K. 

By scrutinizing the level as a function of time I was able to localize 
the effect to the approximate time of the probe removal. Subtraction 
of this effect is covered under data analysis. 

From the above discussion one can see that effects from ground 
radiation are no trivial matter and, in fact, form part of the limit 
on the sensitivity of this whole experiment. This was one of the key 
reasons for observing only the zenith where these effects would be 



minimi zed . 
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B2. THE SWITCH 

In order to study the effect of the switch in detail, I'll first 
consider it as an ideal Faraday rotator. This is done in order to 
determine precisely the signal as a function of the direction of the 



the switch. A convenient reference system for specifying the direction 
of polarization of incoming radiation is the symmetry plane, which 
contains the axis of the ferrite and bisects the right angle formed 
by the orthogonal directions defined by the dual mode transducer (See 
Figure B2.1). Incoming radiation will be specified in terms of the 
magnitude of its electric vector and angle between this vector and the 
symmetry plane 0 o> 

After propagating along the ferrite the E vector has been rotated 
by an amount determined by the longitudinal magnetic field applied to 
the ferrite. In our case we are driving the switch with a sinusoidally 
varying current and thus the angle will be equal to Acoscut + B where 
A is the maximum rotation angle and B is the rotation caused by any 
DC magnetic fields present. Combining these we get an expression for 
the angle between the direction of the E vector and symmetry plane of 
the switch just before the radiation enters the dual mode transducer 



Looking at the signal from one arm of the dual mode transducer, we 
obtain the following for the electric field in the waveguide entering 
the receiver 



incident polarized radiation and of the magnetic field used to drive 



9 = 0 + A coscut + B. 



o 



(B2.1) 




(B2.2) 
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In the detection process the receiver produces a signal proportional 
to the square of the electric field, which, in turn, is proportional 
to the temperature of the received radiation. In this way the signal 
leaving the receiver can be expressed in terms of the temperature of 
the incident radiation as 



Substituting for 9 and expanding, we obtain ' 

CT ( 

S = -g jl + cos2B [sin29 o cos(2A cos(cut)) + cos28 o sin(2A cos(a)t))] 
+ sin2B [cos29 q cos(2A cos(ut)) - sin28 Q sin(2A cos(ut))]j 



for the signal present at the input to the lock-in amplifier , 
where C is a proportionality constant depending on the gain of the 
receiver. 

The action of the lock-in is to extract the Fourier component of 
the fundamental at the switching frequency. The labor can be eased 
considerably by noting that most of the terms in the above expression 
do not have this period. The first term is DC and, of course, gives 
no contribution. The terms in cos (2Acos(aut)) are symmetric about 
at = 0 and have a frequency twice that of the driving current. For 
A = rr/4 the waveform is similar to Figure 2.6b, which also does not have 
a Fourier component at the switching frequency. All that remains are 
terms in sin( 2Acos (at) ), which give signals like Figure 2.6a and are 
the only terms that contribute. By adjusting the phase of the lock-in 
for maximum signal, we account for phase shifts in the apparatus and 
find for the Fourier component and thus the output of the lock-in 




(B2.3) 



(B2.4) 
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' SYMMETRY PLANE 



THE SYMMETRY PLANE 
OF THE SWITCH 



FIG. B2.I 
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V, , . = [cos2B cos20 - sin2B sin20 ] (— ) f 7 sin(2A costut) cosuutdt, 

lock-m 2 o o Joo 4 ' i 

(B2.5) 

where D is a constant depending on the gain of the lock-in. 

62 

The solution to this integral is the Bessel function J^(2A) which gives 0 
V. v • = ^ cos2(9 + B) 2 J {2k) . (B2.6) 

In order to maximize the output, J^(2A) has t^> be maximized. Setting 
dJ^(2A)/dA = 0, we note that the first zero of J^'(x) occurs at x = 1.84l 
and corresponds to a maximum. This implies that for the greatest signal, 
the maximum rotation A is 52.7 degrees. The value for J^(2A) is . 58187 , 
or in terms of the output voltage , 

hock-in = ir (1 - l64) cos2( V B >> < B2 - 7 > 

which differs from that obtained for a purely sinusoidal signal by the 

factor 1.164. This has to be included when interpreting the results. 

The minimum detectable signal for the polarimeter is also altered to 

account for the change in effective amplitude of the modulating signal 
T 

AT . =2/2 — — . (B2.8) 

mn (l.l64)/AvAt 

Returning now to equation B2.6, we note that if the residual fields 
represented by the rotation B are DC in nature, the only effect is a 
change in the effective position of the symmetry plane of the switch, 
an effect which can easily be calibrated out. If there is an AC 
component due to rotation or other movement of the device, we can, 
given that it is small, give the following expression for the lock-in 
output to first order in B , 



120 



V, . . (1.164) [cos26' - 2B' sin20' ], (B2.9) 

lock -in 2 v o o 5 v 

where 6* = 9 + B, and B' is the AC component of B. 
o o dc ^ 

In my case, where the radiometer is mounted vertically and rotated 
about the vertical axis, the maximum possible variation in the 
longitudinal magnetic field is that due to changes in the projection 
of the horizontal component of the earth's field on the axis of the 
ferrite, caused by deviation in the alignments of this axis with the 
vertical. The horizontal component of the earth's field in Princeton 
is about .2 Gauss. Because a misalignment would have a small projection 
in the horizontal direction, .2 Gauss is clearly an upper limit on this 
effect for any orientation. Therefore, an upper limit on the change 
is twice this number. My best estimate of the driving field of the 
switch was 5.3 Gauss and came from wrapping a coil around the switch 
throat and measuring the induced voltage. This number should be less 
than the maximum driving field because the coil was about five centimeters 
above the center of the ferrite, thus it should provide a conservative 
estimate. If one makes the assumption that the Faraday rotation in the 
ferrite is a linear function of field strength, then, the earth's field 
can amount to no more than a 7.5 % effect. In actual practice, magnetic 
shielding of the switch reduced this effect by another factor of 150 
and, hence, to insignificance. 

Another and more important effect and the one which actually made 
the shielding necessary was the change in the absorption coefficient of the 
ferrite with longitudinal magnetization. Figure B2.2a is an abosrption 
curve for a typical ferrite used in microwave work. ° In order to 
demonstrate the effect we are considering, it is only necessary to 



121 



note that, in general, the curve is non-linear and in any particular 
ferrite one is somewhere on the curve in a non-linear region. In a 
device like the switch, where an alternating field is applied, the 
absorption varies sinusoidally with the applied field around the 
zero field value (Figure B2.2b). If we change the zero field 
absorption by biasing the switch with an external magnetic field, 
such as the earth's field, the magnitude of the sinusoidal component 
of the absorption changes due to the change in slope of the curve. 

Because I am detecting the signal with a lock-in amplifier, 

I am only sensitive to changes in the absorption at the switching 
frequency. In order to display this explicitly, I will write the 
absorption, a, in the following form ; 

a = (M) + (M) coscut, (B2.10) 

where the departure from exact sinusoidal dependence to non-linearity 
of the absorption curve has been neglected, as it is only the first 
harmonic I am sensitive to. 

In passing through the switch, unpolarized radiation undergoes 
absorption, which is proportional to the produce, of the absorption 
coefficient and the temperature of the radiation: 



j- = - a T. . (B2.ll) 

ax in v 

In addition, there is a gain in power due to emission in the ferrite, 

giving rise to a signal proportional to the emission coefficient p 

times the temperature of the ferrite T : 




(B2.12) 
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TYPICAL FERRITE ABSORPTION CURVE 





DRIVING MAGNETIZATION 
OF SWITCH 

FIG. B 2.2 (b) 



TYPICAL FERRITE ABSORPTION CURVES 
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The relation between a and j3 can be established from consideration 
of detailed balance. If the input radiation were made to be equal in 
temperature to the switch, the expressions in B2.ll and 12 would have 
to be equal in order that thermal equilibrium be maintained, which 
would imply that O' = p. We can write for the change in temperature 
of the radiation in passing through the switch. 



rim ^ 

E “ ( V V “* < B2 - 13 > 

which has as its solution, 

T .= T. e“ J + T (1 - e -T ) , (B2.l4) 

where t is a times the length of the ferrite and is, in effect, the 
'optical depth' of the switch. For t small this result can be given 
to first order in t by 



T . = (T - T. •) T. 
out o in 



(B2.15) 



Under normal operating conditions T - T. is about 300 K. With the AC 
portion of the absorption this gives an offset signal at the switching 
frequency of 



T = 300 t (M) cosojt°K. (B2.16) 

out/ ac 

Actual measurement of this quantity gave a value of 20 to 100 mdeg. K 

depending on the temperature of the ferrite. This evidence gives two 

pieces of important information. First, it says that t (m) for the 

ac 

-4 

switch is of the order of 3.0 x 10 and thus, direct modulation of the 
input polarized signal by the absorption in the ferrite can be neglected. 
Secondly and more important, it demands that one look carefully at the 
variation of magnetization with radiometer orientation due to the 



\ 
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earth's magnetic field. 

Direct measurement of the sensitivity of the switch to perturbing 
longitudinal magnetic fields showed a change in switch effect of 
156 ± 3 mdeg. K/Gauss for a change of 2.4 Gauss; and a change in 
offset of 130 ±5.5 mdeg. K/Gauss for a change of 1.1 Gauss. 

Application of up to 11 Gauss differences perpendicular to the ferrite 
axis failed to show any effect to the same or^ler of accuracy as the 
longitudinal results. In order to gauge the effect of these offsets, 
assume that the ferrite axis was ten degrees from the vertical. The 
projection of the horizontal component of the earth's field on the 
ferrite would result in a maximum change in the longitudinal field 
of .4 sin (10°) or about .07 Gauss. From the above data this would 
result in a variation in the signal of up to 9.3 mdeg. K, a whopping 
signal! In the early days of debugging the apparatus, signals of this 
order were actually observed and on the basis of this, it was decided 
to shield the switch with high permeability material. Tests run on 
the actual shield used showed that the shield reduced the variation 
in the horizontal component of field by a factor of 150 and all fields 
by better than a factor of 100. In this way the problem of stray magnetic 
fields was also eliminated. 

The final item to be considered is the change in switch offset with 
temperature of the ferrite . Good numbers on this were available from 
the actual data taken in the experiment itself. In addition to measuring 
polarization of the microwave background, I measured the switch 
temperature as a function of time. By averaging together all the data 
at a given temperature, I can obtain the DC offset as a function of 
temperature (see Figure B2.3). It may seem odd that I was able to observe 
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a significant change in temperature in spite of controlling the 
environment inside the experiment shed. As it turned out, heating of 
the horn antenna by the sun during the day and conduction of this heat 
to the switch by the metal waveguide caused an uncontrollable variation 
in the switch temperature. Compensating for extreme changes in 
temperature during a single day was also hard. My controlling apparatus 
had a limited capacity and had to be biased i^ito a proper range by 
controlling the ventilation to compensate for the season of the year. 
Abnormal fluctuation caused the temperature to depart from the controll- 
ing range. 

Luckily, the increase of offset with temperature variation is not 
serious as long as the variations occur on the time scale of a few 
integration times. The analysis process subtracts out offsets and 
under these conditions the temperature induced ones would be eliminated. 
The worst conditions are heating by the sun on a hot day where the 
temperature can go up by as much as 3 degrees c/hr during the heating 
phase and drop by about 2 degrees c/hr in the cooling phase. I 
ignore the heating because data taken during this time is contaminated 
by the sun and not used anyway. The variation during the cooling phase 
and Figure B2.3 can be used to estimate an offset change of approx- 
imately 15 mdeg K/hr during the first couple of hours after the 
sun goes down behind the ground shield. This effect did not 
contribute enough solar signal in the analysis to warrant discarding 



additional data. 



SWITCH OFFSET vs TEMPERATURE 
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FIG. B2.3 
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SPECIFICATIONS OF THE SWITCH 

Ferrotec Rotational Circulator Switch (with 1 dual mode transducer removed) 



Model 


R-165LS 


Serial Number 


11 


Insertion Loss 


.1 db 


Frequency 


8.2-10.0 GHz 


Current for optimum switch rotation 


64 ma 


Optimum switch rotation 


52.7° 


Sensitivity to EL 

Long . 


130 ±5.5 mdeg/Gauss 


Sensitivity to B_l 


0 ± 10.0 mdeg/Gauss 


Offset vs . temp , 


+6.7 mdeg/°C temp. > 3 d°C 




-4.7 mdeg/°C temp. < 36°C 



12’8 



B3. THE RECEIVER 

The receiver consisted of the following components: 
MIXER 

BHG Electronics Laboratory Model MP8/4AARFI 



Serial Number 


5-163-1 


Frequency 


8. 5-9. 6 GHz 


IF Frequency 


60 MHz^ 


IF Bandwidth 


70 MHz 


Power Gain 


29 db 


Noise Figure 
IF AMPLIFIER 


7.2-7. 3 db 



RUG Electronics Laboratory Model EBT101MGC 



Serial Number 


5-163-2 


Center Frequency 


60 MHz 


3 db Bandwidth 


48 MHz 


Maximum Power Gain 
LOCAL OSCILLATOR 


75 db 



Microwave Associates Gunn Oscillator Model MA86101 



Serial Numbex- 


1190 


Frequency 


9.37 GHz 


Power Output 


10 mw 



A useful comparison can be made between the quoted noise 
figure of the mixer given above and the actual system noise. 

A total temperature of l445°K was observed by covering the 
polarimeter with a 294°K absorber and measuring the variance 
of the data. This gives a system temperature of 1445-294 = 1151°K 
It must be noted that this number is referenced to the 290°K 
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temperature difference between the sky and the piece of eccosorb used 
for calibration. Because of insertion loss in the front end of the 
polarimeter, the size of the calibration load is reduced in magnitude 
by the time it reaches the mixer. A rough insertion loss can be 
obtained from summing the following contributions : 



The Switch 
The Circulator 
The Isolator 
Waveguide joints 
Total 



.1 db 
.2 db 
.8 db 

.2 - .3 db 
1.4 db 



Reducing the input reference load by this amount, we find that the 
noise figure is based on a standard of 210 °K rather than 290°K, when 
referenced to the input of the mixer. In these units the system noise 
becomes 

: 833°K. (B3.1) 



T = 1151 
sys 



210 



290 



In addition to absorbing radiation, the components before the mixer 
emit radiation with emission coefficient equal to their absorption 
coefficient. A 1.4 db insertion loss and 290°K device temperatures 
gives rise to 8 l°K worth of emission which must be subtracted from 
the system noise to isolate the mixer contribution: 



T . = 833 - 81 = 753°K. 
mix 

Multiplication by two gives the standard form for a single sideband 



(B3.2) 



T . 
mix 



15o6°k 



= 7.92 db. 



(B3.3) 



This result, when compared with the manufacturer's quoted noise figure 
for the mixer, indicated that there is about a half db worth of extra 
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insertion loss or receiver noise in the system. 

To find the sky noise, a comparison can be made between the 
receiver noise figure and the effective noise figure derived from 
Equation 4.4, using the actual data ffom the experiment. The system 
temperature from the entire years data is about ' l400°K, which is about 
250 degrees higher than the 1151 degree system temperature observed 
with the apparatus terminated in a 300 degree^ load. This additional 
noise is due to the sum of the environmental effects operating on the 
polarimeter and is probably "sky noise" from weather and other 
atmospheric fluctuations. 



b4. tee circulator and isolator 



The factory specifications of the circulator and isolator are: 
THE CIRCULATOR 



Ferrotec Model T-384L 
Serial Number 
Isolation 
Insertion Loss 
VSWR 

Frequency 
THE ISOLATOR 



2 

30.3 db 
.2 db 
1.06 



8.5-10.1 GHz. 



E & M Laboratories Model X112 LI A 



Serial Number 171 

Isolation 40 db 

Insertion Loss .8 db 

VSWR 1.05 



Frequency 



9.2-Q.55 GHz. 
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B5. THE MAGNETIC SHIELDING 

The magnetic shielding was constructed from .050" thick Moly- 
Permalloy, high permeability shielding material by Allegheny Ludlum 
Steel Corporation. It consisted of two 11 inch diameter cylinders each 
20 inches long and so constructed that one would slide within the other. 
Each end was covered with a spun ca.p of the same material, having holes 
for either the antenna and sky horns or the electronics cables. When 
assembled for use, the shield formed a covered cylinder 11 inches in 
diameter and 30 inches long. Moly-Permalloy, when properly annealed 
has the following magnetic properties : 



Maximum Permeability 



Operating Permeability 



85,000-135,000 

300,000-600,000 



Induction at Maximum Perm. 



3000 Gauss 



Saturation Induction 



8000 Gauss 



8.74 



Specific Gravity 
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